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ABSTRACT 



This paper considers a joint replenishment inventory problem with a continuous- 
review (S, c, s) policy for the backorder case with Poisson demands and constant 
procurement lead times. 

Whenever item i's inventory level hits Sj (reorder point) or lower it triggers an 
order so as to raise item i's level to Sj (order up point). At the same time any other 
item j with inventory level at-or-below its can-order point Cj is included in the 
replenishment. 

A Poisson demand model with a queueing description of the system's operation is 
analysed, and comparisions are conducted for joint versus individual orders in the case 
of multi-item problems, where joint replenishment of several items may reduce -setup 
costs. 
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1. INTRODUCTION 



A. DESCRIPTION OF THE PROBLEM 

The purpose in constructing a mathematical model of an inventory system is to 
use it as an aid in developing a suitable operating doctrine for the system. The 
criterion most frequently used for selecting the operating doctrine is that of profit 
maximization or cost minimization. In some cases the task of determining the optimal 
operating doctrine is so difficult that it is either impossible or uneconomical to 
determine the optimal doctrine, and instead, one optimizes with respect to some subset 
of operating doctrines. Occasionally, the mathematical model may be so complicated 
that it is extremely difficult to do anything analytically. In such situations, simulations 
are used to study various operating doctrines. 

We are concerned with finding optimal ordering policies for a multi-item 
inventory’ system in which the inventory position is under continuous review. In this 
thesis, the times at which demands occur are assumed to be generated by independent 
Poisson processes with intensities X- (annual demand rate for item i). We assume that 
all shortages are backordered and will be satisfied from the next shipment. 

We assume that a holding cost of Hj per unit time is charged for each unit of 
item i in inventory. When an order is placed, a setup cost for an independent 

replenishment of item i is assumed to be C- = A + a-, where A is a fixed cost and a^ is 

a cost which depends on the item ordered. The cost of a joint replenishment of both 
item i and j is A + aj + aj . For instance, if item i triggers a replenishment, the setup 
cost of A + aj is applied to item i, and aj to the jointly replenished item j . 

This type of cost structure is particularly appropriate when a group of items is 
ordered from the same supplier, or if a group of items uses the same means of 
transportation, and transportation costs do not increase proportionally when the 

quantity transported is increased. Thus, when a replenishment is made, there is a 

major fixed cost independent of the number of items involved and a variable cost that 
depends upon the number of items [Ref 1: pp. 278-83]. 

Because of this cost structure, it is clear that there may be an opportunity to 
reduce total variable costs if one takes advantage of joint replenishment. 
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B. SCOPE 



Our interest in this problem stems from the "random joint order policy" proposed 
by Balintfy [Ref. 2], and from an (S.c,s) policy treated by Silver [Ref. 3]. For the 
backorder case with Poisson demands and constant lead times, and the specified 
customer service levels, we will compare the individual order policy versus joint order 
policy. 

However, for selecting the control parameters of (S.c,s) in Joint Order Policy, we 
have considered slightly different procedures suggested by Schacck and Silver [Ref 4], 
and Silver [Ref 3]. 

Two programs were written in order to conduct the simulation experiments. The 
first program determines the control parameters of each item in the group for both 
independent and joint replenishment policies. A second program evaluates the model. 
Detailed user instructions and the program itself are contained in .A.ppendix A and B, 
while examples of their use are presented in Chapter II. 

The individual order policy is demonstrated only for the purpose of comparison. 
For the sake of simplicity, we compare the joint order policy to the individual order 
policy where the classical lot size formula is used for each item. 
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II. GENERAL FORMULATION OF THE MODEL 



A. ASSUMPTIONS AND NOTATIONS 

1. Assumptions 

We are dealing with a group of items where the cost of replenishing two or 
more items at the same time is less than the total cost of replenishing the same number 
in separate individual replenishments. The assumptions of the model are: 

1. There is fixed cost A associated uith each replenishment, and variable cost a, 
associated with each item involved in the replenishment. 

2. Demands for item i are Poisson distributed with parameter Xj. 

3. There are no quantity discounts. 

4. Inventorv holding costs are proportional to the average dollar value of 
inventory'. 

5. Service level is defined in one of two wavs: probabilitv of no shortage per 
replenishment cycle, and fraction of demarid to be satis'fied directlv froln the 
shelf. 

2. Basic Notation 

We use the following notation: 

A : fixed cost per replenishment, (independent of the numbei of items and 
units involved) 

a- : fixed cost for including item i in a replenishment. 

Cj : ordering cost per order: ( Cj = A 4- aj) 

EC-: expected relevant annual inventoiy cost for item i. 

Pj : desired customer service level, measured in terms of the probability of no 
shortage per replenishment cycle. 

P 2 desired customer service level, measured in terms of the fraction of 
demand satisfied directly from stock on hand. 

L : replenishment lead time in years, 
r : inventory holding cost rate, 
vj : standard unit price of item i. 

Sj : reorder or must-order point of item i. 
c- : can-order point of item i. 

Sj : order-up-to point of item i. 

Qj : order quantity of item i. 

Xj : Poisson demand rate for item i, in pieces per year. 
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Ij ; average on hand inventory level of item i. 
n : number of items in the replenishment group. 

Nj : expected number of replenishments which involve item i per year. 

NTj : expected number of replenishments triggered by item i per year. 

Other notation will be introduced as they needed. 

B. INVENTORY THEORY BACKGROUND 

Before considering the coordinated control problem it might be helpful to briefly 
discuss the related single item models whose solution will be a key element in 
determing the parameters of the (S, c, s) control system. 

1. Economic Order Quantity Model 

The E.O.Q. model is based on the assumption that the entire lot is added to 
stock at one time, and that the stock will be withdrawn at a constant rate, and no 
stockouts are permitted. 

Total annual variable cost = Ordering cost + Holding cost 



The annual order cost is obtained by multiplying the number of orders per 
year by the cost of placing an order (Cp. 

The average annual holding cost for item i is the average inventory (Qj/2) 
times the annual unit holding cost (Hj= v^r). 

To obtain the minimum cost lot size (E.O.Q.), take the first derivative of total 
annual cost with respect to Q- and set it equal to zero. Solving this for Q^, we get the 
E.O.Q. formula: 



Once the economic order quantity is known, the expected number of orders 
placed during the year, n. can be determined: 



Eq= )^iq/Qi+ HiQi/2 



( 2 . 1 ) 



Qi'-^ = V [2Ciq / Hi] 



( 2 . 2 ) 
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2. Exact Model for the Backorders case >vith Poisson Demands and Constant Lead 
Times 

If demands are random, the probability of stockouts exists, and we must 
include a stockout cost in the expression for expected annual variable costs. This 
expression is given in [Ref. 5: eqn. (4-61)] to be 



where 

: the cost of backorder per unit quantity of backorder. 

Tj : the cost of backorder which is proportional to the length of backorder time. 

Mj : expected lead time demand for item i in pieces, Mj= kjL. 

Hj ; holding cost of item i. 

E(Qi, Sj) : expected number of backorders incurred per year for item i. 

B(Qi, Sj) : the steady state expected number of backorders for item i. 

It is not easy to derive the optimal Qj and Sj by use of the above exact 
expression for ECj. Fortunately, it turns out in practice that it is seldom necessary to 
use the exact formulation, except for the case when it costs very little to incur 
backorders. The approximation discussed in the next section has been shown to give 
good results when backorders are infrequent. 

3. Approximate Form with Known Backorder Cost 



The average annual variable cost equation can be approximated as follows 
[Ref 6: pp. 59 - 62]: 



where 

(Pj : the cost of backorder per unit quantity of backorder. 

E(Mj> Sj) :the expected number of backorders at the end of cycle. 

To minimize ECj, take the partial derivative of the ECj with respect to Qj and 
Sj and set these equal to zero. 




(2.3) 



+ (Tj + Hj)B(Qj, Sj) 




(2.4) 




(2.5) 
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dEC:/ds-= Hj T (<pjXj/ Qj)dE(Mj>Sj)/ds-= 0 



( 2 . 6 ) 



The first equation(2.5) yields the optimal Qj for a given reorder point Sj: 

Qj* = V [2>-i(Ci + (PjE(Mi> S;) / Hj] (2.7) 

The second equation(2.6) yields ; 



dE(Mi>Sj)/dSj = -ZPr(Mj) 



(2.S) 



where 

E (M-> S-) : the expected number of backorders at the end of cycle. 

- ^ ?r (Mj) : complementary cumulative distribution. 

This is the complementary cumulative distribution of the Poisson random 
variable. X, evaluated at Sj. Therefore, solving the second equation for Sj in terms of 
Q- gives the maximum allowable probability of a shortage during a lead time : 

F'(Si) = (2.9) 

To find the optimal pair (Qj, s-) that minimize ECj, iterative procedures can be 
applied, such as described in [Ref. 6: p. 61]. Hower, we do not apply this kind of 
procedure to find out the values of control parameters in (S, c, s) policy, because the 
order size in the (S. c, s) policy is not fixed. 

In (S, c, s) policy, the order size is determined as an order-up-to point minus a 
reorder point for the item which triggers the replenishment action. Order sizes for 
jointly replenishable items are order-up-to points minus on-hand levels which are less 
than or equal to can-buy points but more than reorder points. 

4. Some Useful Formulas for Obtaining the Ii^uted Backorder Cost / Reorder 
Point /Safety Quantity in Fixed Order Size S^>'stem 

In any practical situation, it is very difficult to determine accurately backorder 
costs. They can include such factors as loss of customers' goodwill(i.e., in the future, 
he may take his business elsew'here), or in the military supply system, the cost of 
having some first line weapon system inoperative because of lack of parts. Other 
contributions to shortage cost can be somewhat easier to measure. However, these 
are usually small part of the total backorder costs [Ref. 5: p.l8j. 
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V/e shall assume here that there is a backorder cost associated with each unit 
backordered. 

A complete cost minimization analysis of the optimal values of the model 
parameters can be made only if a cost can be assigned to each shortage. As an 
alternative to specifying the shortage cost, one can examine a given inventory policy 
and obtain the imputed cost of shortage within it [Ref 7: p. 338]. Rather than 
explicitly costing backorders, we will specify a service level constraint and determine 
the imputed stockout cost. 

Whenever an organization uses a service level to establish a reorder point, it 
effectually establishs a stockout cost. Associated with a given service level is an 
imputed or implicit stockout cost. It is a simple matter to determine the imputed 
stockout cost for a given service level from the probability of a stockout [Ref 8; p. 
162). 



a. Probability of no Shortages per Replenishment Cycle 

A service level based on the frequency of service per order interval or 
replenishment cycle will indicate the probability of not running out of stock during the 
replenishment. This approach does not concern itself with how large the shortage is, 
but only with how often it can occur during the lead time (replenishment cycle). It is 
defined as the fraction of the replenishment cycles without depletion of stock: 

number of order periods with a stockout 



Service level fraction per cycle = 1 - 



= l-P,(Mi>Si) 

P,(Mj>Si) = P^(Si) 

= 1 - (service level fraction per cycle) 



total number of order periods 



The term Pj-(sj) is the stockout level fraction per order cycle or the 
probability of at least one stockout while awaiting a supplier's delivery. It is also a 
measure of the fraction of lead time periods during which the demand will exceed the 
reorder point. The magnitude of the stockout is ignored with this approach. 

The following formula is developed to compute the imputed backorder cost: 



Pfs-) = (HjQi) / (bjXj) 



( 2 . 10 ) 



where 

Pj.(s-): probability of stockout for item i. 
Hj : holding cost of item i, (Hj= v^r). 



14 



bj : impuied backorder cost of item i. 



To obtain the appropriate stockout cost, solve eqn.(2. 10) for bj. An 
example will illustrate the procedure in a later section. 

b. Fraction of Demands Satisfied Directly from the Shelf 

The fraction of units demanded and immediately filled from the shelf can be 
defined as the ratio: (number of units supplied)/(total number of units demanded) 
which is equal to (1 - [(number of units short)/(total number of units demanded)]}. 

The above relationships must be measured over some fixed time period 
such as the order interval. 

For example, to obtain the stockout level fraction for units demanded; 

E(M:>S:) 

Stockout level fraction for units = : ^ — 

quantity demanded during a cycle 



E(Mi>Si)/Qi = (HiQi) / (bjAi) 



( 2 . 11 ) 



where 

E(Mj> S-); expected number of backorders for item i during the cycle. 

H- : holding cost of item i. (Hj= v-r). 
bj : imputed backorder cost for item i. 
c. Safety Stock 

With backorders, there is no loss of sales since the customer waits for the 
arrival of the next order if stock is not available. The expected safety stock is defined 
as: 



Safety stock = ^ (s - M) Pj. (M) (2.12) 

= s y Pj. (M) - y M Pj. (M) = s - E(M) 

where 

E(M) ; the expected lead time demand. 

The number of backorders per lead time = 0, if M < s. 

= M - s, if M > s. 
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So, the safety stock is simply the reorder point minus the average lead time 



demand. The reorder point is determined so that the customer service objective is 
satisfied. 

C. MATHEMATICAL STATEMENT OF THE DECISION PROBLEM FOR 
THE JOINT REPLENISHMENT MODEL 

1. Cost Equation 

The basic cost equation in the joint replenishment model is as follow [Ref. 3]: 

EC^= IjVjr + NTjA + (2-13) 

The first term represents the inventory’ holding costs and the other two 
represent the set up charges allocated to item i. The total expected costs per year of the 
group of items are given by 



2. Service Level Constraint 

The service constraint takes one of two forms depending upon the measure of service 
used. 

(1) Probability of no shortage per replenishment cycle 

This is simply the probability that demand for item i during the lead time will 
exceed the reorder point. 



where 

Pj.(s): probability of stock out. 

Mj : lead time demand for item i. 

(2) Fraction of demand satisfied directly from the shelf 

The expected stockout quantity during the lead time is given by 



EC = y ECj= y (I^v-r + NT- A + N-a-) 



(2.14) 



P,(s) = Pr(Mi>Si) = iPr(Mi) 



(2.15) 



E(M > sj) = Z(Mi-Sj)Pr(M-) 



(2.16) 
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where 



E( Mj> Sj) : expected number of stock out for item i in pieces 



during the lead time 
M; : lead time demand for item i 



A backorder occurs if and only if the lead time demand (which is Poisson with 
parameter Mj) is greater than the inventory level of the item when the order is placed. 

D. INDIVIDUAL ORDER POLICY 

Under the the individual order policy, there is a fixed order quantity for item i 
that is ordered everv' time the reorder point is reached . 

Safety stock is needed to protect against stockouts between the time the reorder 
point is reached and the order received. As presented in [Ref. 3], the individual order 
policy (called independent control) is a special case of the joint replenishment policy (S, 
c, s) where c-= sj. In this case item i is ordered only when the inventory' position hits 
its reorder point, and the order quantity is always of size Qj = Sj - Sj. No other items 
will be included in the replenishment, hence a set up cost of A + aj is incurred with 
each replenishment. 

As mentioned in the introduction, the main focus of this thesis is to compare 
joint replenishment with independent control to determine the cost savings that can be 
achieved. 

Let us first consider independent control. Then, Cj = Sj, and the optimization 
problem is: 



Minimize ECj= (Sj - Sj + 1)/ 2 + (Sj - LjL)Vjr 



(2.17) 



+ (A + aj)Lj / (Sj - Sj) 



subject to 



Pj. ( X ^ Sj I XjL) > Pj 



(2.1S) 



or 



1-E(X> Sj|LjL)/Qj> P. 



(2.19) 
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where, X is the lead time demand with a Poisson distribution having mean XjL. 

The left hand side of eqn.(2.lS): Pj.(X ^ Sj| XjL) is the probability of no shortage 
per replenishment cycle. 

In the left hand side of eqn.(2.19), the term E(X>Sj)/Qj is the stockout fraction 
of units demanded. Therefore 1 - E(X>Sj)/Qj is the fraction of demand satisfied 
directly from the shelf 

The exact solution to this problem requires that the pair (Sj, Sj) which minimize 
eqn.(2.17) subject to either eqn.(2.18) or eqn.(2.19) be found for each item i. However, 
it have been observed by many researchers (see: for example Silver and Wilson [Ref 9]) 
that a much simplier sequential approach does almost as well from the standpoint of 
minimizing costs. This sequential approximation approach is what we describe. 

1. Applied Solution Procedure 

Stepl : Let Qj = Sj - Sj. Find the value of .Qj which minimizes ECj ignoring 
the service level constraint. Then given this value of Sj - sj, the constraint is used to 
find the lowest allowable value of Sj. This first step leads to the usual EOQ expression 
Qi = Sj-Sj = V[2(A + aj)Xj Hj], 

Step2 ; Determine the reorder points. 

(1) For the service measure, probability of no shortage per replenishment cycle, 
find the smallest value of Sj, which satisfies the condition of eqn.(2.18): 

Pj.(X ^ SjiXjL) ^ Pj 

(2) For the service measure, fraction of demands satisfied directly from the 
shelf find the smallest value of Sj, which satisfies the following inequality: 

XjL - Sj - XjL Pj. (X < Sj - 1 iXjL) + SjPr (X < Sj I XjL) < Qj(l - P 2 ) 

Steps : Set Sj = Sj + Qj. 

Step4 : Compute the average annual variable costs. 

ECj= (Sj - Sj + 1)/ 2 + (Sj - XjL)vjr + (A + aj)Xj / (Sj - Sj) using the values 
Sj and Sj estimated from Steps 2 and 3. A computer program for solving for these 
parameters is provided in Appendix A. 

Note that we can now obtain imputed stockout costs. 

Steps : Compute the imputed stockout cost 
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(1) For the service measure, probability of no shortage per replenishment cycle: 
Use eqn.(2.10): Pj.(s) = (HjQj) / (bjXj) and solve for b-, bj= H^Qj / (1- P|) Xj. 

(2) For the service measure, fraction of demands satisfied directly from the 
shelf, it is necessary to solve for E(M > Sj) and find the associated Pj.(s) before using 
eqn.(2. 11). It can be done alternatively by finding the complementary cumulative 
distribution of the lead time demand, introduced by eqn. (2.8). Then eqn.(2.11) yields 
bi= HiQi ; [- y (M|)] ,% 

2. Numerical Examples 

Repeated from Silver's example in his paper [Ref. 3]. 

Example 1. Service Level = Probability of no shortage per replenishment 

cycle. 

Consider the following example involving four items: 

A = S 50 

aj = S 10 for i = 1, 2, 3. 4 
r = 0.2 S / S / yr. 

L = I month (i.e. 1/12 year ) 



TABLE I 
INPUT DATA 



item 

(i) 


h 

(piece^'yr) 


(S/pidce) 


P, 

(desired ser\ice level) 


1 


290 


6.90 


0.95 




41 


1.20 


0.95 . 


3 


77 


3.90 


0.95 


4 


122 


2.30 


0.95 



Use of eqn. (2.15), (2.17), (2.18) gives the results shov\m in Table 2. 

Example 2. Service Level = Fraction of demand satisfied directly from the shelf 

Consider the same data as example 1, but with a new the measure of service 
level. Suppose the specified fraction of demands to be satisfied directly is 0.99. The 
results shown in Table 3. 
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TABLE 2 



INDEPENDENT ORDER POLICY: SEPvVICE LEVEL 95% (Pj) 



Item 


LT demand 


EOQ (integer) 


Reorder 


TVC(S) 


S, 0 cost 


I 


24.167 


158.800 (159) 


33 


232.024 


15.113 


2 


3.417 


143.178 (143) 


7 


35.343 


16.762 


3 


6.417 


108.840 (109) 


11 


88.860 


22.051 


4 


10.167 


178.399 a78) 


16 


84.977 


13.453 

1 



TOTAL : 441.203 



TABLE 3 

INDEPENDENT ORDER POLICY: SERVICE LEVEL 99% (P 2 ) 



Item 


LT demand 


EOQ (integer) 


Reorder 


TVC(S) 


S,0 cost 


1 


24.167 


158.800 (159) 


25 


220.984 


*1.983 


2 


3.417 


143.178 143) 


3 


34.383 


1.882 1 




6.417 


108.840 109) 


7 


85.740 


3.499 


4 


10.167 


178.399 178) 


10 


82.217 


1.536 , 



TOTAL : 423.323 



E. JOINT ORDER POLICY 
1. Introduction 

Let us now consider the possibility of joint replenishments. As discussed 
earlier, we can reduce ordering setup costs by combining the orders of several items 
whenever an order is placed. This may result in increased holding costs but the savings 
in ordering costs may be greater than the increase in holding costs resulting in a net 
decrease in total costs. Further, it is likely that the occurance of stockouts will be 
reduced since we will frequently reorders some items before the inventory levels hit the 
normal reorder points for those items. 

We do not seek to determine the "optimal" joint replenishment policy. 
Rather, we restrict attention to a simple, intuitively appealing, policy and seek to find 
the optimal values of parameters for the selected policy. The policy we consider is the 
(S, c, s), "can buy, must buy" policy which requires an order to be placed for item i up 
to the level S- any time the inventory position for item i reaches or falls below s-. In 
addition, however, we order any other item i whose inventory position is below the can 
buy point c-, we always order up to the level S-. 
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2. General Concepts for Obtaining the Values of Control Parameters 

Suppose a single item i, is in a group of coordinated items and assume initially 
that the replenishment lead time is zero. The assumption of unit sized Poisson 
demands together with a zero lead time implies that Sj= 0. 

Now item i, from time to time, will be faced with the opportunity of a 
replenishment at the reduced setup cost of only "aj". Such an event is caused by 
another item triggering a replenishment. Occasionally item i will trigger an order w’hen 
its inventory hits the zero level. Under such circumstances w'e would like to determine 
the (S, c) pair which minimizes the expected costs per unit time. 

Silver shows that the opportunities to replenish at the reduced cost occur 
probabilistically according to a Poisson process with a rate p per year, where ^ is the 
expected number of orders triggered per year by all other items in the group. 

Let NTj be the number of orders triggered per year by item i. Then |tj is given 
by = I NT-, where the summation is taken over all i # j. 

For example, if there are four items in the group and w'e are considering item 
2, we would have =NT^ + NT^ +NT^. 

It is very difficult to determine the value of pj algebraically. We will consider 
some heuristic methods to approximate this. 

Silver [Ref. 3: eqn. 13] determines that the relevant cost equation (the 
subscripts i have been suppressed) is given by: 

Eq= {S - c + p (1 - ) / (1 . p)}-l ((S - c)(S + c + 1) vr / 2 (2.20) 

+ p [c - p (1 - p^ ) / (1 - p)j vr / (1 - p) + X pC A + X a} 
w'here 

p = X / (X + p) 

The algorithm suggested by Silver in [Ref 3] is a simultaneous optimization 
procedure to minimize eqn. (2.17) with respect to Sj, Cj, s-. The procedure is iterative. 

Solution procedures to determine the control parameters of the (S, c, s) policy 
are explained in the next chapter. 
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III. DEVELOPMENT OF THE MODEL 



In this chapter we consider the general joint replenishment problem with positive 
leadtimes. We obtain values of the parameters sj, c-, and Sj using an iterative 
procedure recommended in [Ref. 3]. Silver [Ref 10] derives the following expression for 
the total expected costs per year for item i (we suppress the item subscript for 
simplification of notation): 

Eq = (S-c + p(l -p^)/(l -p)}-l {(S-c)(S + c + 1) vr/ 2 (3.1) 

+ p [c - p (1 - pq/(l - p)jvr/(l . p) + X p^ A + X a} 

. where 

p = X/(X + u) 

p is the probability that a particular occurance is either a demand or an 
opportunity to replenish at reduced cost. 

A. SOLUTION PROCEDURE 

We recommend as a solution procedure a somewhat modified version of the 
procedure suggested by Silver. 

1. Algorithm 

Step 1 : Initialization 

Use the values computed by the independent control to initialize. Let k = 1, 
where k = index of iteration. 

^min = 0- 
^max ~ ^ 

EC(k) = EC(determined from the independent control) 

N(k) = X / q" 

Pj(k) = Nj, where all i # j, and i, j = 1, ..., n 
p(k) = X / (X + p) 
c(k) = 0 
S(k) = q" 
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Step 2 : Let k = k + 1 

c(k) = c(k-l) + / 2, if EC(k) < EC(k-l) 

c(k) = c(k-l) - - W / 2, if EC(k) > EC(k-l) 

S(k) = c(k) - p (1 - p‘^)/(l - p) + V (2 X (a + A X^)/\r 

+ 2cX^'^ ^'(1 - p) - p (1 -p^)(l + p^^*^ ^)/(l - p)2} 

EC(k) = {p (1 - p'^Vd - p)}*^p [c - p (1 - p^'Vd - p)]vr 
/( I - p) + >. p*^ A + X a} 

Step 3 : Test for convergence of c and S. 

Round the c and S values to the nearest integer. 

If c(k) = c(k-l) and S(k) = S(k-l), go to Step 4, else go to Step 2. .And, 
repeat until S and c converge. 

Step 4 ; Repeat Steps 1, 2 and 3 for i = 1, 2, ..., n. 

Step 5 : Let N- = X p ^/[S - c + p (1 - p ^)/(l - p)}. 

p values are from the previous iteration and S and c values are from step 4. 

Step 6 : Repeat steps 1-5 until |c(k) - c(k-l)| <6 and |S(k) - S(k-l)| <5. 

Step 7 : Now determine the must-buy points, Sj, for i = I, ..„ n from the service 

constraints. 

(a) For Service .Vleasure 1: probability (Pj) of no shortage per cycle, find the 
smallest value s such that 

(l/p)Pj.(X<s + c|XL) - Pj.(X=s+l|X L) 

- ps+i y Pj.(X = XglX LHl/pX'^O > 

(b) For Service Measure 2: fraction (P->) of demand to be satisfied directly 
from shelf, find the smallest value s such that 

p*^ {X L - s X L Pj. (X < s - 1 I X L) + s Pj.(X ^ s[X L)} 

+ (1 - P ■'^0 (X L - Wg - X L Pj.(Wg - IIX L) 

+ wp Pj.(X<Wg)} 

< [S-c+pd-p'')/(l-p)}(l-P 2 ) 

Step 8 : Recompute the values Sj and c-. 

Si = S + Si 

C- = C + Sj 

Stop 
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Figure 3.1 Determination of coordinated pair (S,c). 
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IV. SIMULATION STUDY 



In order to determine how costs of ordering, holding, shortage, etc. are affected 
by joint replenishment, a simulation model was developed. The purpose of the 
simulation study was to test the cost savings by comparing the costs for the two 
inventory policies, independent and joint replenishment. The simulation was written in 
FORTRAN??. 

A. ASSUMPTIONS FOR THE SIMULATION STUDY 

1. Poisson demands (exponential interarrival times). 

2. Items demanded one at a time. 

3. Constant lead times (independent of the size of the order and the number of 
items in a replenishment). 

4. For the joint replenishment case, place an order whenever anv item hits its 
must:buv point s- and order up to b-. .Also, include item j # i if the inventorv' 
positionTor item) is below the can-buy point Cj. Order up to Sj. 

5. For the independent control, order up to level S- whenever the stock on hand 

for item i reaches s-. ^ 

6. For a given item, there is never more than a single order outstanding in the 
joint replenishment case. 

B. MODEL INPUTS 

1. Number of items. 

2. Demand rate X- for each item. 

3. Holding cost rate Vj for each item. 

4. Lead time (in months) L. 

5. Type of service level constraint (Typel, or Type2) 

6. Service level required. 

?. Unit price UPj. 

8. Group ordering cost A. 

9. Individual ordering cost aj. 

10. Stockout cost rate per unit stockout (percentage of unit price). 

1 1. Time weighted stockout cost (percentage of unit price). 

12. Sj, Cj, Sj for each item. 

13. Model options (Independent control or Joint Replenishment, etc.) 

14. Limits of simulation period, N'SIM, 
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C. MODEL OUTPLmS 



1. Achieved senice level for each item. 

2. Number of orders (independent and joint). 

3. Total ordering cost for each item. 

4. Total holding cost for each item. 

5. Total unit years of backorders. 

6. Total number of backorders. 

7. Total unit years of stock held. 

S. Total stockout cost (fixed and time weighted) for each item. 

9. Total annual costs for each item. 

10. Total annual variable costs for all items. 

11. Standard deviation of total annual variable costs. 

12. A complete audit trail of all orders placed (optional). 

D. STRUCTURE OF THE PROGRAM 
1. Algorithm 

Step I : Read input data, process options and initialize variables. 

Step 2 ; Determine the earliest event time and type of event. 

Let TDA^ be the time of the next demand for item i and TOAj be the next 
time of the for item of an order for item i. 

Set the master clock time, CT, to the minimum of {TDAj, TOAj}. 

If the event type is a demand for item i then generate another demand time 
for that item. 

If the event type is an order arrival then increament the on hand inventory for 
all items included in the order. 

Step 3 : Determine if an order is required. 

If the onhand level is less than or equal to the reorder point for that item and 
there are no outstanding orders for that item, place an order. If the ordering policy is 
independent, go to step 5 . 

Step 4 : Test for joint replenishments. 

If the onhand level is less than or equal to the can order point for another 
item j in the replenishment set. place joint replenishment for the items to the level Sj. 
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Step 5 : Determine if the simulation is complete. 

If CT ^ NSIM stop, else go to step 2. 

The program collects detailed statistics at each event time for unit years of 
stock held, number of stockouts, unit years of stockouts, number of orders, number of 
units ordered, and number of demands for each item. 

At the conclusion of the simulation additional summar\’ statistics are 
computed for achieved service levels and costs. 
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2. Flowchart 



Init.ia.liza all -raaiabla 



Read input, data & Proceaa Option 




Figure 4.1 Simulation Flowchart. 
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V. EVALUATION OF THE MODEL 

A. PERFORMANCE ANALYSIS 
1. Basic Examples Tested 

Table 4 shows the data used for the evaluation of the joint replenishment 
model. Several combinations (96 cases) of service level, service measure, and lead times 
were considered. The data include the test set considered in Silver's paper [Ref. 3]. 

TABLE 4 

ITEM CHAR.ACTERISTICS 



group 

number 


A 

(S) 


a 

(S) 


number 
of items 


, ITEM DATA 

i A; (piece/yr.) V; (S/piece) 


1 


50 


10 


4 


1 290.000 6.90 

2 41.000 1.20 

3 77.000 3.90 

4 122.000 2.30 


2 


125 


80 


10 


1 40.630 0.45 

2 4.090 36.41 

3 34.680 42.17 

4 4.240 4.46 

5 4.240 4.46 

6 4.090 36.41 

7 4.240 4.46 

8 28.780 3.92 

9 4.090 36.41 

10 4.090 36.41 


3 


50 


5 


4 


1 290.000 6.90 

2 41.000 1.20 

3 77.000 3.90 

4 122.000 2.30 


4 


50 


5 


8 


1 290.000 6.90 

2 41.000 1.20 

3 77.000 3.90 

4 122.000 2.30 

5 50.000 1.20 

6 154.000 3.90 

7 87.000 2.30 

8 25.000 1.20 



2. Model Results 

a. Convergence of Parameter Values and Costs 

Table 5 shows the results of each iteration from Silver's model for joint 
replenishment. Observe that five iterations were required for convergence. This 
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example considers service measure type 2 ■with a specified level 0.99. The lead time is 
one month. 



TABLE 5 

PARAMETERS AND COSTS 



Item 


Iter. 


p 


C 


S 


EC 


1 


1 


0.9942 


79 


146 


201.8498 




2 


0.9972 


84 


152 


210.4767 




j 


0.9973 


84 


152 


210.5555 




4 


0.9974 


87 


152 


210.8846 




5 


0.9973 


87 


152 


210.8239 


2 


I 


0.9272 


31 


90 


21.7035 




2 


0.9499 


36 


97 


23.4369 




j 


0.9464 


36 


96 


23.1011 




4 


0.9469 


36 


96 


23.1429 




5 


0.9467 


36 


96 


23.1271 


3 


1 


0.9650 


27 


85 


66.6963 


• 


2 


0.9768 


48 


89 


70.0377 




3 


0.9741 


43 


88 


68.8164 




4 


0.9744 


43 


88 


68.9584 




5 


0.9744 


43 


88 


68.9309 


4 


1 


0.9774 


45 


138 


63.6090 






0.9850 


71 


145 


66.8895 




j 


0.9835 


70 


143 


65.8678 




4 


0.9835 


70 


143 


65.9083 




5 


0.9835 


70 


143 


65.8689 



b. Cost Savings: Joint Replenishment VS. Independent Control 

Table 6 shows the cost savings achieved by the joint replenishment policy 
for both the modified Silver's method and our heuristic method. 

The comparisions are provided for two cases: 

1. Service measure 1, Pj = 0.95 

2. Service measure 2, P-) = 0.99 

Observe that both the heuristic method and Silver's method achieve cost 
savings greater than 13 % by using joint replenishment, and the heuristic method 
results in greater cost savings than does Silver's method. 

3. Simulation Results 

Figure 5.1 shows the time histoiy’ of stock on hand for items 1 to 4 generated 
by one simulation. The upper graph represents individual control and the lower graph 
represents joint replenishment. 
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TABLE 6 

COMPARISION OF TWO METHODS 



1. modified Silver's method : = 0. 95, L = 1 month 


item 


P 


s 


C 


S 


EC( JP) 


EC( IP) 


Save(%) 


1 


0. 9973 


32 


119 


184 


221. 6339 


232. 024 


4. 48 


2 


0. 9467 


4 


40 


100 


23. 2671 


35. 343 


34. 17 


3 


0. 9744 


9 


52 


97 


70. 9459 


88. 360 


20. 16 


4 


0. 9835 


13 


83 


156 


67. 1722 


84.977 ; 


20. 95 


total cost 1 383.0190 | 441.203 | 13.19 


2. heuristic method : = 0.95, L = 1 month 


item 


P 


s 


C 


S 


EC( JP ) 


EC( IP) 


Save(%) 


1 


0. 9942 


31 


110 


177 


211. 2798 


232. 024 


3. 94 


2 


0. 9272 


3 


34 


93 


21. 6035 


35.343 


3.3. 37 


o 


0. 9650 


9 


36 


94 


63. 7113 


88. 360 


22. 57 


4 


0. 9774 


14 


59 


152 


65. 3723 


84. 977 


23. 07 


total cost 1 351.3467 1 441.203 | 16.83 



3. modified Silver's method : ' 9 ^ = 0.99, L = 1 month 


item 


P 


s 


C 


S 


EC( J?) 


EC( IP) 


Save(%) 


1 


0. 9973 


25 


112 


177 


211. 9739 


220. 984 


4. 08 


2 


0. 9467 


-1 


35 


95 


22. 0671 


34. 383 


35. 82 


3 


0. 9744 


5 


48 


93 


67. 3259 


85. 740 


20. 39 


4 


0. 9835 


7 


77 i 


150 


64. 4122 


82. 217 


21. 66 


total cost 1 366.2788 | 423.323 | 13.48 


4. heuristic method : ?-> = 0.99, L = 1 month 


item 


P 


s 


C 


S 


EC( JP) 


EC( IP) 


Save(%) 


1 


0. 9942 


25 


104 


171 


202. 9998 


220. 984 


8. 14 


2 


0. 9272 


-1 


30 


89 


20. 6435 


34. 383 


39. 96 


3 


0. 9650 


5 


32 


90 


65. 5912 


85. 740 


23. 50 


4 


0.9774 , 


8 


53 , 


146 


62. 6123 


82. 217 


23. 85 


total cost 1 351.8467 | 423.323 | 16.88 
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The simulation results for several combinations of leadtime, service measure 
and service level are represented in Tables 7,8,9. 

When lead time is less than 1 month, the heuristic method yields more savings 
in total cost and also satisfies the required service level. But, as the lead time becomes 
larger, it sometimes fails to meet the required service level, although it does continue to 
generate greater cost savings than Silver's method with the cost parameters examined 
in this thesis. Both methods generate substantial savings over independent control. 
In short lead times, both methods tend to overprotect against stockouts. They 
generally produce service levels in excess of w'hat is required. 

One explanation for the better performance of the joint replenishment models 
for the short lead time examples is the fact that a key step in the joint replenishment 
algorithm for determining the S. c, s parameters assumes that lead time is zero. 
.Adjustments are then made to those values to account for the nonzero lead times but. 
apparently those adjustments are not adequate! .Additional work is required to 
examine this issue and to explore other methods of adjustment. 

We should keep in mind the fact that the model satisfies all requirements in 
the experiment does not mean that it is a correct model — merely that it is a plausible 
one which has not been found inadequate by the data or experiment. 
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ON HAtJD STOCKS ON HAND STOCKS 

BO 120 160 0 40 60 120 160 



2 YEARS RUN : ITEMS IN GROUP 1 





Figure 5.1 Stock on Hand — Independent Control and Joint Replenishment. 



TABLE 7 

COST SAVINGS BY JOINT REPLENISHMENT; WHEN L = 1 MONTH 



1. PARAMETERS FROM MODIFIED SILVER'S METHOD 


Group 

number 


reg. 

SVL 

(%) 


AVG. Total Variable COsts($) 


Cost Savings(%) 
by Joint 
Replenishment 


AVG. SVL 

achieved 

bv 

SIPI. 


INDEPENDENT 


JOINT 


MODEL 


SIM 


MODEL 


SIM 


MODEL 


SIM 


2 


Pl=95 

Pl=99 


441. 20 
448. 06 


443. 04 
445. 77 


383. 02 
392. 74 


361. 48 
368. 70 


13. 19 
12. 35 


IS. 41 
17. 29 


99. 49 
99. 74 


P2=95 

P2=99 


403. 20 
423. 32 


414. 71 
431. 00 


343. 31 
366. 23 


341. 12 ’ 
348. 38 


14. 86 
13. 48 


17. 75 
19. 17 


97. 91 
99. 71 


2 


Pl=95 

Pl=99 


1111. 43 
1153. 33 


1073. 85 
1107. 24 


998. 32 
1037. 70 


953. 78 
956. 40 


10. 18 
10. 03 


11. 22 
13. 62 


99. 27 
100. 00 


P2=95 

P2=99 


1032. 91 

1088. 91 


1029. 70 
1066. 90 


878. 37 
971. 82 


888. 82 
926. 96 


14. 96 
10. 75 


13. 63 
13. 12 


97. 70 
99. 64 


3 


Pl=95 

Pl=99 


423. 30 
430. 16 


422. 05 
425. 74 


356. 89 
365. 33 


345. 93 
354. 15 


15. 69 
• 14. 95 


18. 04 
16. 82 


99. 76 
100. 00 


P2=95 

?2=99 


386. 08 
406. 80 


396. 72 
409. 34 


315. 71 
340. 01 


315. 35 
333. 05 


18. 23 
16. 42 


20. 51 
18. 74 


98. 44 

99. 76 


4 


Pl=95 
Pl=99 1 


676. 10 
686. 94 


677. 55 
683. 89 


513. 14 
528. 48 


490. 92 
508. 801 


24. 10 
23. 07 


27. 54 
25. 60 


99. 21 
99. 80 


P2=9S: 620.56 
P2=99| 650,70 


627. 49 1 
650. 80 1 


458.34 
1 490. 60 


433. 05 26. 14 

459.33 1 24.60 


30. 99 
29. 42 


99. 01 
99. 83 


2. PARAMETERS FROM HEURISTIC METHOD 


1 


Pl=95 

Pl=99 


441. 20 
443. 06 


443. 04 
445. 77 


366. 97 
376. 93 


365. 30 
368. 73 


16. S3 
15. 88 


17. 55 
17. 28 


96. 51 
99. 04 


P2=95 

P2=99 


403. 20 
423. 32 


414. 71 
431. 02 


329. 23 
351. 85 


329. 36 
347. 15 


18. 35 
16. 88 


20. 53 
19. 46 


97. 34 
99. 46 


2 


Pl=95 

?1=99 


1111. 43 
1153. 33 


1073. 85 
1107. 24 


949. 18 
999. 68 


960. 01 
979. 11 


14. 60 
13. 32 


10. 60 
11. 57 


97. 43 
99. 44 


P2=95 

P2=99 


1032. 91 

1088. 91 


1029. 70 
1066. 90 


849. 28 
903. 69 


845. 07 
893. 75 


17. 78 
17. 01 


17. 93 
16. 23 


97. 54 
99. 03 


3 


P"' =Q5 

Pl=99 


423. 30 
430. 16 


422. 05 
425. 74 


332. 99 
343. 17 


3J-2. 7J- 
350. 13 


21. 33 
20. 22 


18. 79 
17. 76 


98. 60 

99. 77 


P2=95 

P2=99 


386. 08 
406. 80 


396. 72 
409. 84 


296. 95 
319. 11 


314. 70 
328. 87 


23. 09 
21. 56 


20. 67 
19. 76 


98. 01 

99. 64 


4 


Pl=95 

Pl=99 


676. 10 
686. 94 


677. 55 
683. 39 


447. 50 
463. 88 


455. 68 
471. 45 


33. 81 
32. 47 


32. 75 
31. 06 


97. 19 
99. 51 


?2=95 

P2=99 


620. 56 
650. 70 


627. 49 
650. 80 


392. 74 
426. 24 


394. 17 
420. 22 


36. 71 
34. 50 


37. 18 
35. 43 


97. 54 
99. 29 



34 



TABLE 8 

COST SAVINGS BY JOINT REPLENISHMENT: WHEN L = 3 MONTH 



1. PARAMETERS FROM MODIFIED SILVER'S METHOD 


Group 

number 


req. 

SVL 

(%) 


AVG. Total Variable COsts($) 


Cost Savings! %) 
by Joint 
Replenishment 


AVG. SVL 

achieved 

bv 

SIM. 


INDEPENDENT 


JOINT 


MODEL 


SIM 


MODEL 


SIM 


MODEL 


SIM 


1 


Pl=95 

Pl=99 


453. 70 
466. 64 


451. 55 
464. 38 


394. 04 
408. 46 


368. 34 
383. 34 


13. 15 
12. 47 


18. 43 
17. 45 


98. 97 

99. 74 


P2=95 

P2=99 


404. 20 
431. 48 


407. 47 
433. 39 


340. 27 
373. 06 


331. 60 
350. 99 


15. 82 
13. 54 


18. 62 
19. 01 


96. 91 
99. 59 


2 


Pl=95 

?1=99 


1173. 94 
1224. 36 


1135. 01 
1177. 53 


1019. 79 
1078. 73 


938. 40 
994. 54 


13. 13 
11. 89 


17. 32 
15. 54 


98. 33 

99. 33 


P2=95 

P2=99 


1047. 78 
1109. 34 


1054. 40 
1097. 29 


893. 13 
992. 50 


855. 98 
917. 02 


14. 76 
10. 53 


18. 82 
16. 42 


98. 34 

99. 78 


3 


Pl=95 

Pl=99 


435. 80 
, .448. ,74 


439. 10 
446. 89 


365. 75 
382. 33 


358. 36 
371. 04 


16. 07 
. 14. 80 


18. 39 
16. 97 


96. 73 
99. 51 


P2=95 

P2=99 


386. 76 
413. 82 


390. 60 
410. 99 


314. 71 
347. 03 


309. 32 
337. 40 


18. 63 
16. 14 


20. 81 
17. 91 


97. 88 
99. 57 


4 


Pl=95 

iPl=99 


1 694. 66 
713. 84 


690. 92 
711. 89 


525. 31 
549. 89 


491. 81 
516. 17 


24. 38 
22. 97 


23. 82 
27.49 1 


99. 32 
99. 80 


P2=95 

?2=99 


619. 70 
660. 82 


623. 99 
660. 38 


446. 77 
498. 85 


418. 44 
467. 03 


27. 91 
24. 51 


32.94 1 

29.28 


98. 70 

99. 84 


2. PARAMETERS FROM HEURISTIC METHOD 


1 


Pl=95 

Pl=99 


453. 70 
466. 64 


451. 55 
464. 38 


377. 98 
392. 64 


363. 73 
378. 46 


16. 69 
15. 86 


19. 45 
18. 50 


97. 35 
99. 76 


P2=95 

P2=99 


404. 20 
431. 48 


407. 47 
433. 39 


326. 78 
357. 32 


318. 15 
346. 98 


19. 15 
17. 19 


21. 92 
19. 94 


96. 75 
99. 44 


2 


Pl=95 

Pl=99 


1173. 94 
1224. 36 


1135. 01 
1177. 53 


978. 83 
1041. 41 


966. 47 
1015. 42 


16. 62 
14. 94 


14. 85 
13. 77 


*93. 78 
99. 17 


P2=95 

P2=99 


1047. 78 
1109. 34 


1054. 40 
1097. 29 


854. 13 
950. 04 


360. 07 
923. 59 


13. 48 

14. 36 


18. 43 
15. 83 


97. 94 
99. 44 


3 


Pl=95 

Pl=99 


435. 30 
448. 74 


439. 10 
446. 89 


340. 91 
358. 43 


344. 09 
360. 88 


21. 77 
20. 13 


21. 64 
19. 25 


99. 20 
99. 77 


P2=95 

P2=99 


386. 76 
413. 82 


390. 60 
410. 99 


291. 57 
324. 41 


301. 49 
327. 00 


24. 61 
21. 61 


22. 32 
20. 44 


97. 70 
99. 74 


4 


Pl=95 

Pl=99 


694. 66 
713. 84 


690. 92 
711. 89 


456. 20 
483. 56 


447. 29 
474. 04 


34. 33 
32. 26 


35. 26 
33. 41 


95. 34 
98. 97 


P2=95 

P2=99 


619. 70 
660. 82 


623. 99 
660. 38 


381. 62 
432. 52 , 


400. 09 
425. 76 


38. 42 
34. 55 


35. 88 
35. 53 


95. 55 
99. 10 
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TABLE 9 

COST SAVINGS BY JOINT REPLENISHMENT: WHEN L = 6 MONTH 



1. PARAMETERS' FROM MODIFIED SILVER'S METHOD 


Group 

number 


req. 

SVU 

(%) 


AVG. Total Variable COsts($) 


Cost Savings! %) 
by Joint 
Repleni shment 


AVG. SVL 

ac^Hved 

bv 

SIM. 


INDEPENDENT 


JOINT 


MODEL 


SIM 


MODEL 


SIM 


MODEL 


SIM 


1 


Pl=95 

Pl=99 


465. 47 
485. 50 


461. 92 
482. 24 


404. 32 
425. 84 


374. 87 
396. 24 


13. 14 
12. 29 


18. 85 
17. 83 


97. 53 
99. 24 


P2=95 

?2=99 


406. 25 
441. 41 


407. 04 
440. 34 


343. 22 
381. 50 


320. 37 
354. 81 


15. 51 
13. 57 


21. 29 
19. 42 


96. 23 
99. 34 


2 


Pl=95 

Pl=99 


1217. 73 
1285. 89 


1174. 03 
1240. 05 


1063. 46 
1131. 83 


983. 25 
1045. 79 


12. 26 
11. 98 


16. 25 
15. 67 


97. 90 
99. 33 


P2=95 

P2=99 


1046. 00 
1161. 47 


1056. 59 
1132. 14 


897. 43 
1014. 36 


851. 39 
933. 91 


14. 20 
12. 67 


19. 42 
17. 51 


98. 07 

99. 77 


3 


Pl=95 

Pl=99 


447. 56 
467.. 60 


442. 91 
464. 06 


375. 58 
397. 10 


359. 91 
. 380. 81 


16. 08 
15. 08 


18. 74 
17. 94 


94. 30 
99. 14 


?2=95 

P2=99 


390. 96 
423. 50 


392. 25 
418. 99 


317. 44 
355. 48 


210. 39 
342. 46 


18. 81 
16. 06 


20. 87 
18. 27 


97. 19 
99. 29 


4 


Pl=95 

|P1=99 


712. 21 
740. 67 


711. 50 
737. 77 


536. 53 
571. 58 


505. 25 
529. 95 ■ 


24. 66 
j 22.83 


28. 99 
28. 17 


96. 47 
99. 04 




P2=95 

|P2=99 


624. 71 
674. 31 


629. 13 
674. 72 


447. 04 
507. 72 


411.24 28.44 

465.69 1 24.76 


34.64 j 97.19 
30.98 1 99.44 


2. PARAMETERS FROM HEURISTIC METHOD 


1 


Pl=95 

Pl=99 


465. 47 
485. 50 


461. 92 
482. 24 


387. 43 
409. 79 


379. 63 
401. 59 


16. 77 
15. 60 


17. 82 
16. 72 


*92. 54 
98. 74 


P2=95 

P2=99 


406. 25 
441. 41 


407. 04 
440. 34 


328. 41 
366. 23 


328. 97 
363. 02 


19. 16 
17. 03 


19. 18 
17. 56 


95. 33 
99. 08 


2 


Pl=95 

Pl=99 


1217. 73 
1285. 39 


1174. 03 
1240. 05 


1030. 96 

1090. 96 


991. 45 
1049. 95 


15. 34 
15. 15 


15. 55 
15. 33 


98. 92 
100. 00 


P2=95 

P2=99 


1046. 00 
1161. 47 


1056. 59 
1132. 14 


848. 43 
962. 32 


878. 96 
940. 00 


18. 89 
17. 15 


16. 81 
16. 97 


96. 53 
99. 34 


3 


Pl=95 

Pl=99 


447. 56 
467. 60 


442. 91 
464. 06 


351. 19 
372. 95 


341. 63 
362. 02 


21. 53 
20. 24 


22. 37 
21. 99 


-91. 90 
*96. 14 


P2=95 

P2=99 


390. 96 
423. 50 


392. 25 
418. 99 


292. 41 
331. 47 


298. 23 
330. 47 


25. 21 
21. 73 


23. 97 
21. 13 


95. 57 
98. 48 


4 


Pl=95 

?1=99 


712. 21 
740. 67 


711. 50 
737. 77 


466. 69 
503. 17 


532. 09 
559. 74 


34. 47 
32. 07 


25. 21 
24. 13 


*36. 62 
*94. 34 


P2=95 

P2=99 


624. 71 
: 674.81 


629. 13 
674. 72 


381. 93 
442. 41 


404. 55 
452. 60 


38. 36 
34. 44 


35. 70 
32. 92 


*92.36 1 

*97. 70 



36 



VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

For each of several test combinations, we determined the independent and 
coordinated control parameter values and the associated expected costs. We simulated 
several years of operations with the derived joint replenishment ordering parameters to 
evaluate the effectiveness of the models as compared to independent control models. 
We also use the simulation to compare a simple heuristic joint replenishment model. 
With the simulation we were able to accumulate costs (ordering, holding and stockout) 
and effectiveness as measured by two service levels. 

In summarv’, the following statements can be made: 

1. Coordinated control results in total costs sisnificantlv lower than under 
independent control. The average cost savings are about fS.T % . 

2. Coordinated control provides better service than independent control. In a few 
cases, the actual service for an individual item is slishtlv lower than required; 
when this occurs, the item involved is the one that tnsgers the majority of the 
replenishments in the group. 

3. Cost savings of coordinated control increases as the ratio a/A decreas es. 

4. Cost savincs imorove as the number of items in replenishment set increases; the 
more items there are in the group, the more beneficial coordinated 
replenishment becomes. 

5. For a given service measure and lead time, the cost savings increase as the 
desired service level decreases. 



B. RECOMMENDATIONS 

The procedure to determine the parameters for coordinated control requires only 
the data for independent control. We observe that the joint replenishment models 
perform best for shorter lead times. This is probably the result of the fact that the 
joint replenishment algorithms initially assume zero lead times and then adjust the 
derived parameters heuristically to account for nonzero lead times. This adjustment 
procedure may be inadequate. Additional w'ork is needed to consider other adjustment 
methods. 
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APPENDIX A 

"PARA": PRO GRAM TO SELECT PARAMETERS 



VARIABLE DEFINITION ** 



DM 

MR 

S/Sl 

C/Cl 

IPC 

IPS 

tci 

EC/ECl 

TLM 

TLT 

TL 

DLT 

NIT 

TN 

R / RHO 
F 

UP 

OCI 

OCG 

SVL 

SCU 

PSO 

STON 



ANNUAL DEMAND 
REORDER POINT 
ORDER - UP POINT 
CAN ORDER POINT 
OPTIMAL PAIR OF C 
OPTIMAL PAIR OF S 
ECNOMIC ORDER QUANTITY 

AVERAGE ANNUAL TOTAL VARIABLE COST (INDEPENDENT ODER) 
AVERAGE ANNUAL TOTAL VARIABLE COST(JOINT ORDER) 

LEAD TIME (IN MONTH) 

LEAD TIME (IN YEARS) 

AVERAGE LEAD TIME (IN YEAR) 

LEAD TIME DEMAND 

NO. OF TOTAL ITEM IN THE GROUP 

NO. OF TRIGGERING AN ORDER BY ITEM I IN THE GROUP 
OPPORTUNITY TO REPLENISH AT REDUCED COST 
HOLDING COST RATE TO PROCUREMENT COST 
UNIT PRICE 

INDIVIDUAL ORDERING COST 
GROUP ORDERING COST 
SERVICE LEVEL 

IMPUTED UNIT STOCKOUT COST GIVEN SVL. 

MAX. ALLOWABLE ROBABILITY OF STOCK OUT 
EXPECTED NUMBER OF STOCKOUT DURING A ORDER CYCLE 



DIMENSION ITM(IOO) ,DM(100) ,UP(100) ,H( 100 ) , OCI ( 100 ) ,ECI(100) ,Q(100) 
&,C(100) ,TN(100),S(100) ,IS(100) ,R(100) ,U(100) ,EC(100) ,MR(100) , 
&IQ(100) ,DLT(100) ,SCU(100) ,IC(100) ,RHO(50,100) ,EC1(50,100) , 
&C1(50,100) ,31(50,100) ,IC1 (50, 100) ,IS1( 50, 100), IPC( 100) , IPS (100) , 
&MR1(100) ,TL(100) 

* DATA INITIALIZATION 
STN=0 . 

TVC=0 . 

* READ PROCESS OPTIONS 

CALL OPTION ( ID , LT , KEY , LTT , SVL , MSRCH ) 

READ(ID,110)NIT,OCG,F 

TLT=LT/12. 

DO 11 1=1, NIT 

READ(ID,111)ITM(I) ,DM(I) ,UP(I) ,OCI(I) 

^ H(I)=UP(I)*F 
IF(LTT.EQ.1)TL(I)=TLT 
IF(LTT.EQ.2)THEN 
CALL LEAD(TLT,TLY) 

TL(I)=TLY 
END IF 

11 CONTINUE 

* PRINT SELECTED OPTIONS 

WRITE ( 7 , 1 00 ) KEY , SVL , LTT , ID , MSRCH 

* CALCULATE E . 0 . Q 

DO 1 I =1,NIT 
DLT(I)=DM(I)*TL(I) 
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OC=OCG+OCI(I) 

0(I)=SgPsT(2*DM(I)*0C/H(I)) 
fQ I =INT Q(I +.5) 



* COMPUTE REORDER POINT :IND. ORDER POLICY 



KEY = 1 : PROS. OF NO SHORTAGE PER REPLENISHMENT CYCLE 

KEY =2': FRACTION OF DEMAND TO BE SATISFIED DIRECTLY FROM SHELF 



CALL POISON(KEY,DLT(I) ,0(1 ) ,SVL,MX,M,MAX,CCDF1) 

MR1(I)=M 
IF(KEY.EQ.1)THEN 

SCU( I )=H(I)*Q(I )/((!.- SVL)*DM(I)) 

END IF 

IF(KEY.EQ.2)THEN 

SCU(I)=H(I)*Q(I)/(CCDF1*DM(D) 

END IF 
ECI(I)=H(I) 

WRITE(7,501 
TVC=TVC+ECI 
TN(I)=DM(I)/Q(I) 

STN=STN+TN(I) 

1 CONTINUE 

WRITE(7,502)TVC 

* DETERMINE THE PARAMETERS FOR JOINT ORDER POLICY 

(REORDER / CAN ORDER / ORDER-UP POINT) 

DO 2 1=1, NIT 

2 U(I)=STN-TN(I) 

DO 3 1=1, NIT 

3 R(I)=DM(I)/(DM(I)+U(D) 

* SEARCH LOOP FOR DETERMINIG A PAIR (S,C) : ITEM BY ITEM 



*( (Q(I)+1.)/2.+(MR1(I)-DLT(I)))+0C*(DM(I)/Q(I)) 
)lTM(lhDLT(I),TL(I),Q(I),IQ(I),MRl(I),ECI(l),SCU(I) 
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5 



11=1 

SEC=0. 

SEC1=0. 

DO 4 1=1, NIT 

RH0(I,II)=R(I) 

K=1 

CMIN=0 . 

CMAX=Q(I) 

C(K)=0. 

S(K)=Q(I) 

EC(K)=ECI(I) 

TEST=10.**10 

ICASE=1 

LOOP=0 

K=K+1 

IF(ICASE.EQ.1)THEN 
C(K)=C(K-l)+(CMAX-CMIN)/2 
END IF 

IF(ICASE.EQ.2)THEN 
C(K)=C(K-lh(CMAX-CMIN)/2 
END IF 

S(K) = C(K)-R(I)*(l-R(I)**C(Kn/ 
& +OCG*R(I)^*C(K))/H 



(1-R(I))+SQRT(2*DM(I)*(0CI(I) 

..a) + 2*C(K)’^R(I)**(C(K) + l)/(l-R(I) 
R(I)* 1-R(I)*^C(K))*(1+R I ** C(K)+l| j/(l-R I))**2 



EC(K)=(1./(S(K)-C(K)+R(I)*(1.-R(I)**C(K|)/(1.-R(I)))) 
&*((S(K)-C(K))*(S(K)+C(K)+l)*H(I)/2+R(I)^(cU)-Ru)^(l •-R(I)**C(K) 
&/(l.-R(l)))=^H(I)/(l.-R(I))+DMa)*R(I)*^C(X)*OCG+DM(I)’*'OCI(I)) 



) 



;k)=int(s 

k)=int(c 



(K)+isj 



is 

IC 

IF(C(K) .EO.C(K-l) .AND.S(K) .EQ.S(K-1))G0 TO 6 

IF(EC(K' 

IF(EC( 



K). 

K). 



GT. (0.9995^EC(K-D) .AND.EC(K) .LE.EC(K-D) GO TO 6 
LT.EC(K-l)) THEN 
IF(C(K) .GT.C(K-1))THEN 
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7 

4 

20 

21 

22 

C 

23 



CHIN=C(K) 

ICASE=: 

GO TO 5 
ELSE 

CMAX=C(K) 
ICASE=2 
GO TO 5 
END IF 



IF(C(K) .LE.C(K-1))THEN 
CMIN=C(K) 

ICASE=1 
GO TO 5 
ELSE 

CMAX=C(K) 

ICASE=2 
GO TO 5 
END IF 

END IF 

IF(LOOP.EQ.O)THEN 

CMIN=0 

CMAX=C(2) 

ICASE=2 
L00P='l 
GO TO 5 
END IF 

* FIND OUT THE PAIR (S,C), MINIMIZES EC 
DO 7 J=1,K 

IF(EC(J) .LE.TEST)THEN 
TEST=EC(J) 



CONTINUE 

ICl(I,II)=INT(Cl(I,II)+.5) 

IS1(I ,Il)=INT(Sl(I,II)+.5) 

SEC1=SEC1+EC1(I,II) 

CONTINUE 

SHAVE=(TVC-SEC1)/TVC 

WRITE(7, 507)11 
DO 20 1=1, NIT 

WRITE(7,508)I,RHO(I,II) ,IC1(I,II) ,IS1(I,II) ,EC1(I,II) 

CONTINUE 

WRITE (7, 509 )SEC1, SHAVE 

IF(MSRCH.EQ.2) GO TO 42 

* COMPUTATIONS NEEDED ONLY FOR S. METHOD 1 

IF(II.GE.2)THEN 
DO 21 1=1, NIT 

IF(C1(I,II).EQ.C1(I,II-1).AND.S1(I,II).EQ.S1(I,II-1)) GO TO 42 
END IF 

IF(II.GE.2)THEN 
DO 22 1=1, NIT 

IF(II.GT.NIT.0R.EC1(I,II) .GT. (0.9995*EC1( 1,11-1 )) 

& .AND.EC1(I,II) .LE.EC1(I,II-1)) GO TO 42 

END IF 

IF (II.GE.7) GO TO 42 
DO 23 1=1, NIT 

TN(I) = DM(I) * R(I)**C1(I,II) / ( Sld.Il) - C1(I,II) + R(I) 

& *(1 - R(I)-^Cl(I,Ilh / (1 -R(l))) 

CONTINUE 
SUM=0 . 



ELSE 



EC1(I,II)=EC(J) 




END IF 
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24 

25 

26 



C 

42 



27 



C 

C 

C 



28 



C 

100 



110 

111 

200 

501 

502 

503 

504 

505 



DO 24 1=1 .HIT 
SU1I=SUM+TN(I) 

DO 25 :=1.MIT 
ua)=SUH-TN(l) 

DC 23 1=1, HIT 
R(I)=DH(I)/(DM(I)+U(I)) 
11 = 11+1 
SEC1=0. 

GO TO 41 



WRITE(7,503) 

DO 27 1=1, NIT 
C(lj=ClfI,II) 
S(I)=S1(I,II) 
IPC(Ij=INT(C(I)+.5) 



IPS(l)=INT(S(l)+.5) 
WRITE(7, 504)1, IPC(I), 
WRITE (7, 200) 



IPS(I) 



* COMPUTE REORDER POINT: JOINT ORDER POLICY 



DO 28 1=1, NIT 
R(I)=RHO(I,II) 

IF(KEY.EQ.1)THEN 

CALL MUST1(R(I) ,C(I) ,DLT(I) ,SVL,MBUY) 
MR(I)=MBUY 
C(I)=C(I)+MR(I) 

S(I)=S(I)+MR(I) 

IC(I)=INT(C(I)+.5) 

IS(I =INT(S(l)+.5) 

EC(I)=EC1(I,II)+(MR(I)-DLT(I))*H(I) 

END IF 

IF(KEY.E0.2)THEN 

CALL flUST2(R(I) ,C(I) ,S(I) ,DLT(I) ,SVL,MBUY) 
MR(I)=MBUY 
C(I)=C(I)+MR(I) 

S(l)=SaKMR(I) ^ 

IC(I =INT(C(I)+.5) 

IS(l)=INT(S(I)+.5) 

EC(I)=EC1(I,II)+(MR(I)-DLT(I))*H(I) 

END IF 

WRITE(7,505)I,R(I) ,MR(I) ,IC(I) ,IS(I) ,EC(I) 
SEC=SEC+EC(I) 

SHAVE= (TVC-SEC) /TVC 
WRITE (7 , 506 ) SEC , SHAVE 
STOP 



FORMAT ( IX, /, '** SELECTED OPTIONS : ' , / , IX, ' (* KEY = ',12 , 2X, ' * SVL. 
& =' ,F4.2,2X, '* LEAD TIME TYPE =' ,I2,2X, '* DATA=' ,12, 

&2X,'* S.METHD=' ,12, IX, ')',//, ' * KEY= 1 : SVL. BY PROB. OF NO SHOR 
STAGE PER REPLENISHMENT CYCLE',/,' 2 ; SVL. BY FRACT. OF DEM 

TO BE SATISFIED DIRECTLY FROM SHELF ',/,' * LTT= 1 : CONSTANT 
&LEADTIME 2 : RANDOM (UNIFORM W/ MEAN) * DATA=1 : DATA 

&1 2 : DATA2 3 : DATA3 4 : DATA4 ',/,' * SEARCH METHOD =1 : 

& MODIFIED SILVER METHOD 2 : HEURISTIC, 

&///,20X, '*** INDEPENDENT ORDER POLICY ***' 

&,//,lX, 'ITEM' ,3X, 'LT DEMAND(LT.) ' ,2X, 'EOQ(N-INT. ) ' ,2X, ' (INTEGER) ' , 
&2X, 'REORDER' ,2X, 'EXP. TVC($) ' ,2X, ' S/0 COST') 

FORMAT ( I 3 , 2X , F6 . 2 , 2X , F4 . 2 ) 

F0RMAT(I3,2X,F7.3,3X,F5.2,2X,F5.2) 

FORMAT(//,20X, '*** JOINT ORDER POLICY ***',/, 2X, ' ITEM' , 5X, ' RHO ' 
& , 8X , ' REORDER ' , 5X , ' CANBUY ' , 5X , ' ORDERUP ' , 5X , ' EXPECTED COST ' ) 
F0RMAT(1X,I3,4X,F7.3, ' ( ' , F5 .3 , ' ) ' , 3X, F8 . 3 , 6X, 14 , 7X, 13 ,4X, F9 .3 , 3X, 
83) 

FORMAT (3 IX, 'SUM OF AVG. ANNUAL TVC($)= ' , F9 . 3) 

FORMAT( ' ' ,// ,15X, OPTIMAL PAIR (C,S) ',/, 

&16X, 'ITEM' ,7X, 'OPT. C',7X,'OPT. S') 

FORMATa6X,I3,9X,I4,9X,I4) 

FORMAT(2X,I3,4X,F6.4,7X,I4,9X,I4,7X,I4,7X,F10.4) 
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506 FORMAT U5X, 'TOTAL COST =' ,F10.4, ' ( ' ,F6.4, ' ) • ) 

507 FORtlAT( " ,// ,3X. '* PARAMETERS & EK. COSTS FOR JOINT REPLENISHMENT- 
&,/,15X, ' (V^TiEN LEADTIME = 0)’,/,l5X,'* AT ITERATION :',I2,/, 

&1X, 'ITEM' ,4X, 'RHO' ,4X, 'OPT.C ,2X, 'OPT. S' ,3X, 'EXP. COST' ) 

508 FORMAT(2X,I2,3X,F6.4,4X,I3,4X,I3,4X,F8.4) 

509 FOPilAT(23X, 'TOTAL =' ,F9.4, ' ( ' ,F6.4, ' ) ' ) 

END 

C 

SUBROUTINE POISON (KEY,DLT,Q,SVL,MX,M,MAX,CCDF1) 

C 

PSO=l-SVL 
IF(KEY.EQ.l) THEN 

CALL SERV1(DLT,PS0,MAX,M) 

RETURN 
END IF 

IF (KEY. EQ. 2) THEN 
CALL SERV2(Q,DLT,PS0,MAX,M,CCDF1) 

RETURN 
END IF 
END 
C 

SUBROUTINE MAXD(DLT,MAX) 

C 

DIMENSION PMF(750) 

1=1 

PMF(1)=EXP(-DLT) 

CDF=PMF(1) 

GO TO 2 

1 PMF(I) = (DLT/(I-1))'^PMF(I-1) 

CDF=CDF+PMF(I) 

2 IF(CDF.GE. 0.9995) GO TO 3 
1 = 1+1 

GO TO 1 

3 MAX =1+1 
RETURN 
END 

C 

SUBROUTINE SERVl (DLT , PSO ,MAX ,M) 

C 

DIMENSION PMF(200) ,CDF(200) ,CCDF(200) 

PMF(1)=EXP(-DLT) 

CDF(1)=PMF(1) 

CCDF(1)=1-CDF(1) 

CALL MAXD(DLT,MAX) 

DO 1 I=2,MAX+1 

PMF(I) = (DLT/(I-1) )=''PMF(I-1) 

CDF(I)=CDF(I-1)+PMF(I) 

1 CCDF(I)=1-CDF{I) 

DO 2 I=1,MAX+1 

2 IF(PSO.LT.CCDF(I) .AND.PSO.GT.CCDF(I+l))GO TO 3 

3 M=I 
RETURN 
END 

C 

SUBROUTINE SERV2 (Q , DLT , PSO , MAX ,M , CCDFl ) 

C 

DIMENSION PMF(200) ,CDF(200) ,CCDF(200) 

STON=Q*PSO 

PMF(1)=EXP(-DLT) 

CDF(1)=PMF(1) 

CCDF(1)=1.-CDF(1) 

CALL MAXD(DLT,MAX) 

DO 1 I=2,MAX+1 

PMF(I)=(DLT/(I-1))*PMF(I-1) 

CDF(I)=CDF(I-1)+PMF(I) 

CCDF(I)=1-CDF(I) 

1 CONTINUE 
1=0 

2 IF(I.EQ.O)THEN 
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PST=DLT 

"“fST=DLT- I - DLT*CDF ( I ) + 1 *CDF ( I t 1 ) 

END IF 

IF(PST.LT.STON) THEN 
K=I 

CCDF1=CCDF(M+1) 

ELSE 
1 = 1+1 
GO TO 2 
END IF 
RETURN 
END 
C 

SUBROUTINE LEAD(TLT,TLY) 

C 

C 

CALL LRND(719325,U, 1,1,0) 

TLY=2*TLT*U 

RETURN 

END 

C 

SUBROUTINE MUSTl ( R , C , DLT , S VL , MBUY ) 

C 

DIMENSION PMF(400) ,CDF(400) 

IC=INT(C+.5) 

RC=(1./R)’^’^IC 

RHS=SVL/(R**IC) 

PMF(1)=EXP(-DLT) 

IF(PMF(1) .LT.l.OE-70) PMF(1)=0. 
CDF(1)=PMF(1) 

DO 1 1=2,400 

PMF(I)=(DLT/(I-1))*PMF(I-1) 

IF(PMF(I) .LT.l.OE-70) PMF(I)=0. 
CDF(I)=CDF(I-1)+PMF(I) 

1 CONTINUE 
JS=-2 

3 SUM=0 . 

KF=JS+2 
KL=JS+IC 
DO 2 IX=KF,KL 

2 SUM=SUM+PMF(IX+1)*((1/R)**IK) 

TEST=RC*CDF ( KL+1 ) -PMF ( JS+2 ) - (R** ( JS+1 ) ) *SUM 
IF(TEST.GE.RHS) THEN 

MBUY=JS 

ELSE 

JS=JS+1 
GO TO 3 
END IF 
RETURN 
END 
C 

SUBROUTINE MUST2 (R , C , S , DLT , SVL , MBUY) 

C 

DIMENSION PMF(400) ,CDF(400) 

IC=INT(C+.5) 

RHS=(S-C+R*(1-R**C)/(1-R))*(1-SVL) 

PMF(1)=EXP(-DLT) 

IF(PMF(1) .LT.l.OE-70) PMF(1)=0. 
CDF(1)=PMF(1) 

DO 1 1=2,400 

PMF(I)=(DLT/(I-1) )*PMF(I-1) 

IF(PMF(I) .LT.l.OE-70) PMF(I)=0. 
CDF(I)=CDF(I-1)+PMFU) 

1 CONTINUE 

JS=-1 

3 SUM=0 . 

DO 2 IW=JS+1, JS+IC 
IF(IW.LT.l) THEN 
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VCDF=0 . 

IF( IW.EQ.O)VCDrl=CDF(l) 

ELSE 

VCDF=CDF(IW) 

VCDF1=CDF(IW+1) 

END IF 

2 SUM = SUM+(1/R**IW)*(DLT-IW-DLT*VCDF+IW*VCDF1) 

IF(JS.LT.l) THEN 
VCDF=0 . 

IF( JS.EQ.0)VCDF1=CDF(1) 

ELSE 

VCDF=CDF(JS) 

VCDF1=CDF(JS+1) 

END IF 

TEST=R**C*(DLT-JS-DLT*VCDF+JS*VCDF1)+(1-R)*R**(C+JS)*SUM 
IF(TEST.LE.RHS) THEN 
MBUY=JS 
ELSE 

JS=JS+1 
GO TO 3 
END IF 
RETURN 
END 
C 

SUBROUTINE OPTION (ID,LT,KEY,LTT,SVL,MSRCH) 

C 

WRITE (6, 50) 

READ(5,51) ID 
WRITE (6, 53) 

READ(5,51) key 
WRITE(6,55) 

READ (5, 57) SVL 
WRITE(6,52) 

READ(5,51) LT 
WRITE (5, 54) 

READ (5, 51) LTT 
WRITE(6,56) 

RE AD ('5,51) MSRCH 
RETURN ' 

C 

50 FORMATCl',' * ENTER OPTION FOR DATA FILE ( 1- 4) ! 1 = 

&DATA1 2 = DATA2 3 =DATA3 4 = DATA4') 

51 FORMAT(Il) 

52 FORMATS 1',' * ENTER THE LEADTIME IN MONTH ! ',/,5K, '1 = 1 MONTH 

&',/,5X,'3 = 3 MONTH' ,/ ,5X, '6 = 6 MONTH') 

53 FORMATCl',' * ENTER OPTION FOR THE METHOD OF SVL. MEASURE ! ',/,' 

& 1 = PROB. OF NO SHORTAGE PER REPLENISHMENT CYCLE ',/,4X,' 2 

&= FRACTION OF DEMAND TO BE SATISFIED DIRECTLY FROM SHELF ' ) 

54 FORMATCl',' ENTER THE LEADTIME TYPE ! ',/,5X, '1 = CONSTANT 
& ',/,5X,'2 = RANDOM (UNIFORM W/ MEAN)') 

55 FORMATCl',' * ENTER THE DESIRED SERVICE LEVEL ! ',/ ,4X, ' (SHOULD 
&BE IN THE FORM OF "F4.2" : EX. ; 0.90 OR 0.99 ...ETC.)') 

56 FORMATCl',' * ENTER OPTION FOR SEARCH METHOD U-2) !',/,' 1 = 

&MODIFIED SILVER METHOD 2 = HEURISTIC) 

57 FORMAT (F4. 2) 

END 



1. INPUT FORMAT 



004 


050.00 


0.20 




1 


290.000 


6.90 


10.00 


2 


41.000 


1.20 


10.00 


3 


77.000 


3.90 


10.00 


4 


122.000 


2.30 


10.00 
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2. SXSC PROGRAM 



&TRACE OFF 
ScFN = PARA 

&FNO = &CONCAT OF &FM OUTPUT 

&TYPE Do you need to compile your program ? (Y) 

&READ VAR &R_COMPILE 

&IF &R_COMPILE NE Y &GOTO -RUN 

-H FORTVS &FN 

&IF &RC EQ 0 &GOTO -RUN 

ScTYPE Your program did not compile; check for errors. 

&TYPE Do you wish to view the program LISTING file? (Y) 
&READ VAR &RSP1 

&IF &RSP1 EQ Y BROWSE &FN LISTING A 

ScTYPE Do you wish to XEDIT the program file? (Y) 

ScREAD VAR ScRESPl 

ScIF ScRESPl NE Y ScEXIT 1 

ScCOMMAND XEDIT ScFN FORTRAN A 

ScTYPE Do you wish to run the program again? (Y) 

ScREAD VAR ScRESP2 

ScIF ScRESP2 EQ Y ScGOTO -H 

ScEXIT 1 

-RUN 

FILEDEF 01 DISK ScFN DATAl A1 

FILEDEF 02 DISK ScFN DATA2 A1 

FILEDEF 03 DISK ScFN DATA3 A1 

FILEDEF 04 DISK ScFN DATA4 A1 

FILEDEF 07 DISK ScFN OUTPUT A1 (LRECL 133 

LOAD ScFN (START 

ScIF ScRC EQ 0 ScSKIP 9 

ScTYPE Your program did not run correctly; check for errors. 
ScTYPE Do you wish to XEDIT the program file? (Y) 

ScREAD VAR ScRESP3 

ScIF ScRES?3 NE Y ScEXIT 2 

ScCOMMAND XEDIT ScFN FORTRAN A 

ScTYPE Do you wish to run the program again? (Y) 

ScREAD VAR &RESP4 

ScIF ScRESP4 EQ Y ScGOTO -H 

ScEXIT 2 

ScTYPE YOUR OUTPUT IS IN THE FILE ScFN OUTPUT A 
ScTYPE Do you wish to BROWSE your output? (Y) 

ScREAD VAR ScRESP 

ScIF ScRESP EQ Y &COMMAND BROWSE ScFN OUTPUT A 
ScTYPE Print your output file? (Y) 

&READ VAR ScRESP? 

ScIF ScRESP? EQ Y ScCOMMAND PRINT ScFN OUTPUT A 
-REDO 

ScTYPE Do you wish to XEDIT the program file? (Y/N) 

ScREAD VAR &RESP5 

ScIF &RESP5 EQ Y XEDIT ScFN FORTRAN A 

ScTYPE Do you wish to run the program again? (Y) 

ScREAD VAR &RESP6 

ScRESPSe = ScCONCAT OF ScRESPS ScRESPS 
ScIF ScRESPSe EQ YY ScGOTO -H 
ScIF ScRESPS EQ Y ScGOTO -RUN 
S^IXIT 
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3. RUN EXAMPLE 



S ara 

o you need to compile your program ? (Y) 

?S FORTRAN COMPILER ENTERED. 11:58:02 

**MAIN** END OF COMPILATION 1 ****** 

**POISON** END OF COMPILATION 2 ****** 

**MAXD** END OF COMPILATION 3 ****** 

**SERV1** END OF COMPILATION 4 ****** 

**SERV2** END OF COMPILATION 5 ****** 

**LEAD** END OF COMPILATION 6 ****** 

**MUST1** END OF COMPILATION 7 ****** 

**MUST2** END OF COMPILATION 8 ****** 

**OPTION** END OF COMPILATION 9 ****** 

VS FORTRAN COMPILER EXITED. 11:58:10 



EXECUTION BEGINS . . . 

* ENTER OPTION FOR DATA FILE ( 1- 4) ! 

1 = DATAl 2 = DATA2 3 =DATA3 4 = DATA4 

1 



* ENTER THE LEADTIME IN MONTH ! 
1=1 MONTH 
3=3 MONTH 
6=6 MONTH 



* ENTER OPTION FOR THE METHOD OF SVL. MEASURE ! 

1 = PROB. OF NO SHORTAGE PER REPLENISHMENT CYCLE 

2 = FRACTION OF DEMAND TO BE SATISFIED DIRECTLY FROM SHELF 

2 



* ENTER THE LEADTIME TYPE ! 

1 = CONSTANT 

2 = RANDOM (UNIFORM W/ MEAN) 

1 



* ENTER THE DESIRED SERVICE LEVEL ! 

(SHOULD BE IN THE FORM OF "F4.2" : EX. ; 0.90 OR 0.99 ...ETC.) 
0.99 

* ENTER OPTION FOR SEARCH METHOD (1-2)! 

1 = MODIFIED SILVER METHOD 2 = HEURISTIC 

1 

YOUR OUTPUT IS IN THE FILE PARA OUTPUT A 
Do you wish to BROWSE your output? (Y) 

^rint your output file? (Y) 
n 

Do you wish to XEDIT the program file? (Y/N) 
n 

Do you wish to run the program again? (Y) 
n 

R; T=1.47/2.39 11:58:39 
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4. OUTPUT 1 ( P2 = 0.99 ) 



SELECTED OPTIONS : 

(* KEY = 2 * SVL. =0.99 * LEAD TIME TYPE = 1 * DATA= 1 * S.METHD= 1 ) 

* KEY= 1 : SVL. BY PROB. OF NO SHORTAGE PER REPLENISHMENT CYCLE 

2 : SVL. BY FRACT. OF DEMAND TO BE SATISFIED DIRECTLY FROM SHELF 

* LTT= 1 : CONSTANT LEADTIME 2 ; RANDOM (UNIFORM W/ MEAN) 

* DATA=1 : DATAl 2 : DATA2 3 : DATA3 4 : DATA4 

* SEARCH METHOD =1 : MODIFIED SILVER METHOD 2 : HEURISTIC 



ITEM 


*** INDEPENDENT 
LT DEMAND (LT.) EOQ(N-INT.) 


ORDER POLICY *** 

(INTEGER) REORDER 


EXP. TVC($) 


S/0 COST 


1 


24.167(0.083 


1 158.800 


159 


25 


220.984 


1.983 


2 


3.417(0.083 


1 143.178 


143 


3 


34.383 


1.882 


3 


6.417(0.083 


1 108.340 


109 


7 


35.740 


3.499 


4 


10.167(0.083! 


1 178.399 


178 


10 


82.217 


1 . 536 






SUM OF 


AVG. ANNUAL 


TVC($) 


= 423.323 





* 


PARAMETERS 


& EX 
(WHEN 
^ AT 


.COSTS FOR 
LEADTIME ^ 
ITERATION 


JOINT REPLENISHMENT 
= 0) 

: 1 


ITEM 


RHO 


OPT.i 


C OPT.S 


EXP. COST 


1 


0.9942 


79 


146 


201.3498 


2 


0.9272 


31 


90 


21.7035 


3 


0.9650 


27 


85 


66.6963 


4 


0.9774 


45 


138 

TOTAL = 


63.6090 

353.3584(0.1641) 



* PARAMETERS & EX. COSTS FOR JOINT REPLENISHMENT 







(WHEN 


LEADTIME 


= 0) 






^ AT 


ITERATION 


: 2 


ITEM 


RHO 


OPT.i 


C OPT.S 


EXP. COST 


1 


0.9972 


84 


152 


210.4767 


2 


0.9499 


36 


97 


23.4369 


3 


0.9768 


48 


89 


70.0377 


4 


0.9850 


71 


145 


66.3895 








TOTAL = 


370.8406(0.1240) 



* PARAMETERS & EX. COSTS FOR JOINT REPLENISHMENT 
(WHEN LEADTIME = 0) 

^ AT ITERATION : 3 



ITEM 


RHO 


OPT.C 


OPT.S 


EXP. COST 


1 


0.9973 


84 


152 


210.5555 


2 


0.9464 


36 


96 


23.1011 


3 


0.9741 


43 


88 


68.8164 


4 


0.9835 


70 


143 

TOTAL = 


65.8678 

368.3406(0.1299) 


* ] 


PARAMETERS 


& EX. COSTS FOR 
(WHEN LEADTIME = 
* AT ITERATION 


JOINT REPLENISHMENT 
= 0) 

: 4 


ITEM 


RHO 


OPT.C 


OPT.S 


EXP. COST 


1 


0.9974 


37 


152 


210.8846 


2 


0.9469 


36 


96 


23.1429 


3 


0.9744 


43 


88 


68.9584 


4 


0.9835 


70 


143 

TOTAL = 


65.9083 

368.3940(0.1286) 



47 



* PARAMETERS & EX. COSTS FOR JOINT REPLENISHMENT 
(tffiEN LEADTIME = 0) 

^ AT ITERATION s 5 



ITEM 


PwHO 


CPT.C 


OPT.S 


EXP. COST 


1 


0.9973 


67 


152 


210.8239 


2 


0.9467 


36 


96 


23.1271 


3 


0.9744 


43 


88 


68.9309 


4 


0.9835 


70 


143 

TOTAL 


65.8689 

= 368.7505(0.1289) 



** OPTIMAL 


PAIR (C,S) 


** 


ITEM 


OPT. C 


OPT. : 


1 


37 


152 


2 


36 


96 


3 


43 


88 


4 


70 


143 







*** JOINT 


ORDER POLICY 


•k-kiK 




ITEM 


RHO 


REORDER 


CANBUY 


ORDERUP 


EXPECTED COST 


1 


0.9973 


25 


112 


111 


211.9739 


2 


0.9467 


-1 


35 


95 


22.0671 


3 


0.9744 


5 


48 


93 


67.8259 


4 


0.9835 


7 


77 


150 

TOTAL COST 


64.4122 

= 366.2788(0.1348) 
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5. OUTPUT 2 (PI = 0.95 ) 



** SELECTED OPTIONS : 

(* KEY = 1 * SVL. =0.95 * LEAD TIME TYPE = 1 * DATA= 1 * S.METHD= 1 ) 

* KEY= 1 : SVL. BY PROB. OF NO SHORTAGE PER REPLENISHMENT CYCLE 

2 : SVL. BY FRACT. OF DEMAND TO BE SATISFIED DIRECTLY FROM SHELF 

* LTT= 1 ; CONSTANT LEADTIME 2 : RANDOM (UNIFORM W/ MEAN) 

* DATA=1 : DATAl 2 : DATA2 3 : DATA3 4 : DATA4 

* SEARCH METHOD =1 : MODIFIED SILVER METHOD 2 : HEURISTIC 



ITEM 


*** INDEPENDENT 
LT DEMAND (LT.) EOQ(N-INT.) 


ORDER POLICY *** 

(INTEGER) REORDER 


EXP. TVC($) 


S/0 COST 


1 


24.167(0.083 


) 158.800 


159 


33 


232.024 


15.113 


2 


3.417(0.083 


> 143.178 


143 


7 


35.343 


16.762 


3 


6.417(0.083 


» 108.840 


109 


11 


88.860 


22.051 


4 


10.167(0.083; 


) 178.399 


178 


16 


84.977 


13.453 






SUM OF 


AVG. ANNUAL 


TVC($) 


= 441.203 





* 


PARAMETERS 


& EX. COSTS FOR 
(WHEN LEADTIME ^ 
* AT ITERATION 


JOINT REPLENISHMENT 
= 0) 

; 1 


ITEM 


RHO 


OPT.i 


C OPT.S 


EXP. COST 


1 


0.9942 


79 


146 


201.3498 


2 


0.9272 


31 


90 


21.7035 


3 


0.9650 


27 


85 


66 . 6963 


4 


0.9774 


45 


138 

TOTAL = 


63.6090 

353.3584(0.1980) 


* 


PARAMETERS 


& EX 
(WHEN 
^ AT 


.COSTS FOR 
LEADTIME • 
ITERATION 


JOINT REPLENISHMENT 
= 0) 

: 2 


ITEM 


RHO 


OPT.i 


C OPT.S 


EXP. COST 


1 


0.9972 


84 


152 


210.4767 


2 


0.9499 


36 


97 


23.4369 


3 


0.9768 


48 


89 


70.0377 


4 


0.9850 


71 


145 

TOTAL = 


66.8895 

370.8406(0.1595) 



* PARAMETERS & EX. COSTS FOR JOINT REPLENISHMENT 







(WHEN 
^ AT : 


LEADTIME 

ITERATION 


= 0) 
: 3 


ITEM 


RHO 


OPT.I 


C OPT.S 


EXP. COST 


1 


0.9973 


84 


152 


210.5555 


2 


0.9464 


36 


96 


23.1011 


3 


0.9741 


43 


88 


68.8164 


4 


0.9835 


70 


143 

TOTAL = 


65.8678 

368.3406(0.1651) 



* PARAMETERS & EX. COSTS FOR JOINT REPLENISHMENT 







(WHEN LEADTIME 
^ AT ITERATION 


= 0) 
: 4 


ITEM 


RHO 


OPT.C 


OPT.S 


EXP . COST 


1 


0.9974 


87 


152 


210.8846 


2 


0.9469 


36 


96 


23.1429 


3 


0.9744 


43 


88 


68.9584 


4 


0.9835 


70 


143 

TOTAL = 


65.9083 

368.8940(0 
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* PARAMETERS & EX. COSTS FOR JOINT REPLENISHMENT 
(WHEN LEADTIME = 0) 

* AT ITERATION : 5 



ITEM 


RHO 


OPT.C 


OPT.S 


EXP. COST 


1 


0.9973 


87 


152 


210.8239 


2 


0.9467 


36 


96 


23.1271 


3 


0.9744 


43 


88 


68.9309 


4 


0.9835 


70 


143 

TOTAL 


65.8689 
= 368.7505 



** OPTIMAL PAIR (C,S) ** 

ITEM OPT. C OPT. S 



1 


37 


152 


2 


36 


96 


3 


43 


88 


4 


70 


143 







*** JOINT 


ORDER POLICY 


'k’k:k 




ITEM 


RHO 


REORDER 


CANBUY 


ORDERUP 


EXPECTED COST 


1 


0.9973 


32 


119 


184 


221.6339 


2 


0.9467 


4 


40 


100 


23.2671 


3 


0.9744 


9 


52 


97 


70.9459 


4 


0.9835 


13 


83 


156 


hi .1122 



TOTAL COST = 383.0190(0.1319) 
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oonnnnnnooonnnnnnooooooonnoonnooonnnnnonnoonnonnonno 



APPENDIX B 

"SIM'':PROGRAM TO TEST THE MODEL 



VARIABLE DEFINITIONS 



KEY = OPTIONS FOR ORDER POLICY 

ISTA = OPTIONS FOR START INVENTORY LEVEL 

IPRINT= OPTIONS FOR PRINTOUT 

NIT = NO. OF ITEMS IN THE GROUP TO BE EVALUATED 

NSIM = SIMULATION PERIOD 

SIMT = SIMULATION INTERVAL TIME (lYEAR) 



DM = ANNUAL DEMAND RATE 
ME = ORDER-UP POINT 

lEQ = ECONOMIC ORDER QUANTITY( INTEGER) 

MC = CAN ORDER POINT (CAN BUY) 

MR = REORDER POINT(Ml=IND . :M2=JOINT) 

TDA = TIME OF DEMAND OCCURING TIME 
TOA = ORDER QTY. ARRIVAL TIME 
M = INDEX OF EVENT FOR DEMAND OCCURED ITEM 

lA = INDEX OF EVENT FOR OREDR QTY. ARRIVED ITEM 

lET = INDEX OF EVENT TYPE 

CT = MASTER CLOCK TIME IN YRS(THE EAliLIEST EVENT OCCURING TIME) 

DA = TIME INTERVAL BETWEEN 2 SUCCESSIVE DEMANDS 

ARV = UNIFORM RANDOM VARIABLE FOR LEADTIME WITH MEAN "ULT" 

ULT = EXPECTED LEADTIME 
IQ = ORDER QTY. 

IT = INDEX FOR REPLENISH TRIGGER ITEM 

JT = INDEX FOR JOINT REPLENISHED ITEM 

lOH = ONHAND INVENTORY 

lOHS = START LEVEL 

UOH = UNIT YEARS ONHAND QTY 

UBO = UNIT YEARS BACKORDER QTY 

NBQ = NO. OF BACKORDERS (QTY.) 

NTBQ = TOTAL NUMBER OF BACKORDERS 

BOT = BACKORDER TIME 

TBOT = TOTAL BACKORDER TIME 

NOD = TOATL NUMBER OF ORDERS FOR EACH ITEM 

NDM = TOTAL NUMBER OF DEMANDS FOR EACH ITEM 

NTDM = TOTAL NUMBER OF DEMANDS DURING TMAX 

NJTO = NUMBER OF JOINT REPLENISHMENTS 



UP = UNIT PRICE 

F = HOLDING COST RATE 

HC = HOLDING COST 

OC = ORDER & REVIEW COST 

OCI = INDIVIDUAL ORDERING COST 

OCG = JOINT (GROUP) ORDERING COST 

SCU = BACKORDER COST PER BACKORDER UNIT 

SCT = BACKORDER COST (TIME DEPENDENT FACTOR) 

TVC = AVG. ANNUAL TOTAL VARIABLE COSTS 



DIMENSION ISEED(IO) ,DM(10) ,MR(10) ,MC(10) ,ME(10) ,IOH(10) ,ASVL(10), 
&TDA(100) ,IT(10,100), JT(10,100),IQ(100) ,DA(100) ,TOA(100) ,IEQ(10), 
&H(10) ,OCI(10) ,UPao) ,SCU(10) ,SCTa0),F(10),NJTOh00) ,OCT(100) , 
&UOH(10,100) ,AUOH(lOO) ,UBO(10,100) , AUBO(100 ) ,NOD( 10 , 100 ) ,KB ( 10 , 100 ) 
&,ANOD(100),NDM(10,100) ,ANDM(100) ,NBQ(10,100) ,NTBQ(10 , 100) ,KBT(10) , 
&BOT(10 ,100) ,TBOT(10,100) ,FRC(10,100) , FRO( 10 , 100 ) , SVL( 10 , 100 ) , 
5AFRQ( 100) ,00(10,100) ,AOC(100) , COST( 10 , lO'O ) ,TCOST ( 100) ,NTOD(100) , 
&HCOST(10,100) ,SCOST(10,100),SCTU(10,100),SCTT(10,100) ,NTDM(100) , 
&AHCOST(100) ,ASCTU(100) ,ASCTT(100) ,AB0T(10) ,ABOQ(100) ,AC0ST(100) , 



51 



&raCOST(100) ,TSCTU(100) ,TSCTT(iOO) ,OCTT(100) ,K1(10) ,M2(10) , 
&TUCH(100KTUBO(100 LtFRT(10) ,TFRQ(100),NTDMS(100) ,NTCDS(100) , 
SAN3(10) 

COMMON I CASE , I S EED , DM , DA , ARV , ULT , I , I STP , MR , ME , KEY , lEQ 
C 

C DATA INITIALIZATION 

C ■ ■ • 

SIMT = 1. 

TMAX = 0 . 

IC =0 
ISD = 293715 
CT =0. 

DT =0. 

AT =0. 

K =1 

N =0 

TVC = 0 . 

NTJTO = 0 
TANDM = 0 . 

TANOD = 0 . 

TAOC = 0 . 

TCST = 0. 

STV = 0 . 

ANJ = 0. 

C 

C READ PROCESS OPTIONS 

C 

CALL OPTI ON ( ID , RATE 1 , RATE 2 , KEY , LTOPT , NS IM , I STA , IP ) 

C 

C READ INPUT DATA 

C 

READ(ID,101)NIT,OCG,ISV,MSV,LT 

ULT=LT/12. 

DO 401 1=1, NIT 

READ(ID,100)DM(I) ,M1(I) ,M2(I) ,MC(I) ,ME(I) ,IEQ(I) ,OCI(I) ,UP(I) ,F(I) 
IF (KEY . EQ . 1 )MR ( I ) =M1 ( I ) 

IF(KEY.E0.2)MR(I)=M2(I) 

SCU(I)=RATE1*UP(I) 

SCT(I)=RATE2*UP(I) 

H(I)=UP(I)*F(I) 

401 CONTINUE 
C 

WRITE (7 , 198 ) ID , KEY , NSIM , ISTA , IP , LTOPT , MSV , ISV , ULT , LT 
C 

DO 1 1=1, NIT 

ISEED(I)=ISD+I**3 
KBT(I) = 0 
ANB(I) = 0. 

TFRT(I) = 0. 

ABOT(I) = 0. 

ASVL(I) = 0. 

TOA(I) = 10. **10 
1 TDA(I) = 0. 

DO 2 1=1, NIT 
DO 2 J=l,100 

DA(J) = 0. 

IQ(J) = 0 

SVL(I,J) = 0. 

UOH(I,J) = 0. 

TUOH(J) = 0. 

AUOH(J) = 0. 

COST(I,J) = 0. 

TCOST(J) = 0. 

ACOST(J) = 0. 

UB0(I,J) =0. 

TUBO(J) = 0, 

AUBO(J) = 0. 

NOD(I,J) = 0 

NTOD(J) = 0 
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NTODS(J) 

ANOD^J) 

NDMd, J) 

NTDH(J) 

HTDMS(J) 

ANDH(J) 




0 

0 . 

0 

0 

0 

0 . 

0 

0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 , 

0 . 

0 . 

0 . 



OCTT(J) = 0. 

AOC(J) = 0. 

NJTO(J) = 0 

HCOST(I,J) = 0. 

THCOST(J) = 0. 

AHCOST(J) = 0. 

SCOST(I,J) = 0. 

SCTU(I,J) = 0. 

TSCTU(J) = 0. 

ASCTU(J) = 0. 

SCTT(i,J) = 0. 

TSCTT(J) = 0. 

ASCTT(J) = 0. 

KB (I,J) = 0 

IT (I,J) = 0 

JT (I,J) = 0 

DO 3 1=1, NIT 

IF(ISTA.EQ.1.AND.KEY.EQ.1)THEN 
lOH(I) = lEQ(I) 

END IF 

if(ista.eq.i.and.key.eq.2) then 

lOH(I) = ME(I) 

END IF 

IF(ISTA.EQ.2) THEN 
CALL START 
IOH(I)=ISTP 
END IF 
CONTINUE 

PRINTOUT INPUT DATA 

DO 4 1=1, NIT 

WRITE(7,199) I,DM(I),IEQ(I) ,MR(I) ,MC(I) ,ME(I) ,IOH(I) ,ULT 
WRITE (7, 97) 

DO 5 1=1, NIT 

WRITE (7, 98) I,UP(I),F(I) ,H(I) ,OCG,OCI(I) ,SCU(I) ,SCT(I) 

IC = IC + 1 

IF (IC.GE. 65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2) THEN 
IC = 0 

WRITE(7,200) 

END IF 
CONTINUE 

-- START SIMULATION 



IF (IP.EQ.1.0R.IP.EQ.2)THEN 
WRITE (7, 200) 

END IF 



DO 777 ITER=1,NSIM 

IF (ITER.GT.1,AND.IP.EQ.1.0R.ITER.GT.1.AND.IP.EQ.2) THEN 
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IC = 0 

WRITE n, 200) 

END IF 

TMAX=T1IAX+SIMT 

20 IF(CT.GT.TMAX) GO TO 77 
N = N+1 
JOINT = 0 

DETERMINE THE INITIAL DEMAND OCCURANCE TIMES FOR EACH ITEM 

IF(N.LE.l) THEN 
DO 6 1=1, NIT 

ICASE=I 
CALL DEMAND 
TDA(I)=DA(ICASE) 

6 UOH(I,ITER)=UOH(I,ITER)+IOH(I)*DA(ICASE) 

DETERMINE THE EARLIEST EVENT 

(DEMMAND OCCURED / ORDER PLACED / ORDERED QTY. ARRIVED) 

IF(TDA(1) .LE.TDA(2)) THEN 
TD=TDA(1) 

ELSE 

TD=TDA(2) 

END IF 

DO 7 1=1, NIT 
IF(TDA(I) .LE.TD) THEN 

td=tdaU) 

M=I 

END IF 

7 CONTINUE 

IET=0 

CT=TDA(M) 

IOH(M)=IOH(M)-l 

IF(MSV.E0.1.AND.I0H(M) .EQ.-l) KB(M,ITER)=KB(M, ITER)+1 
NDM(M,ITER)=NDH(M,ITER)+1 

SET NEXT DEMAND OCCURING TIME FOR THE ITEM ;'M' 

ICASE=M 

CALL DEMAND 

TDA(M)=CT+DA(ICASE) 

GO TO 90 
END IF 
C 

IF(TDA(1) .LE.TDA(2)) THEN 
TD=TDA(1) 

ELSE 

TD=TDA(2) 

END IF 

DO 8 1=1, NIT 

IF(TDA(I) .LE.TD) THEN 
TD=TDA(I) 

M=I 
END IF 

8 CONTINUE 
C 

IF(TOA(l) .LE.TOA(2)) THEN 
TA=TOA(l) 

ELSE 

TA=TOA(2) 

END IF 

DO 9 1=1, NIT 

IF(TOAd) .LE.TA) THEN 
TA=TOA(I) 

IA=I 
END IF 
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9 CONTINUE 
C 

IF(TA.LE.TD)THEN 
GO TO 93 
ELSE 
IET=0 
CT=TD 
ICASE=M 
CALL DEMAND 
TDA(M)=CT+DA(ICASE) 

IOH(M)=IOH(M)-l 

IF(MSV.E0.1.AND.I0H(M) .EQ.-l) KB (M , ITER)=KB(M, ITER)+1 
NDM(M,ITER)=NDM(M,ITER)+1 
END IF 



DETERMINE SHOULD ORDER 



90 



C 



IF(IOH(M) .GT.MR(M))THEN 

UOH(M,ITER)=UOH(M,ITER)+IOH(M)*DA(M) 

IF(IP.EQ.1)THEN 

WRITE (7, 501 )CT,IET,M,UOH(M, ITER) , (lOH(I) , 1=1 ,4) 

END IF 

IF(IP.EQ.2)THEN 

WRITE (7, 601 )CT,IET,M,UOH(M, ITER) ,IOH(M) 

END IF 

IC = IC + 1 

. IF (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.6-5.AND.IP.EQ.2) THEN 

IC = 0 

WRITE(7,200) 

END IF 
GO TO 20 
END IF 

IF(K.EQ.1)THEN 

IET=1 

K=K+1 

IT(M,K)=M 

IF(IOH(M) .LT.O) THEN 
NBQ(M,ITER)=-IOH(M) 

IF(KEY.EQ.l) IQ(M) = IEQ(M)+NBQ(M, ITER) 

IF(KEY.EQ.2) IQ(M) = ME(M)+NBQ(M,ITER) 
IF(LTOPT.EQ.l)TOA(M) = CT + ULT 
IF(LTOPT.EQ.2)THEN 
CALL ARRIVE 
TOA(M)=CT+ARV 
END IF 

BOT(M,ITER)=TOA(M) - CT 

IF(NBQ(M,ITER) . EO . 1 )TBOT(M , ITER)=T30T(M, ITER)+BOT(M, ITER) 
UBO (M , ITER)=UBO (II , ITER)+BOT (M , ITER ) 

IF(IP.EQ.1)THEN 

WRITE(7, 502)CT,IET,M,NBQ(M, ITER) ,BOT(M, ITER) ,UBO(M,ITER) , 
& IQ(M) ,TOA(M) , (I0H{lhl=l,4) 

END IF 

IF(IP.EQ.2)THEN 

WRITE(7,602)CT,IET,M,NBQ(M,ITER) ,BOT(M, ITER) , UBO(M, ITER) , 
& IQ(M) ,TOA(M),IOH(M) 

END IF 

IC = IC + 1 

IF (IC.GE.65.AND.IP.EQ.1 .OR.IC.GE.65.AND.IP.EQ.2) THEN 
IC = 0 

WRITE(7,200) 

END IF 
ELSE 

IF(KEY.EQ.l) IQ(M) = lEQ(M) 

IF(KEY.EQ.2) IQ (M)=ME (M) -IOH(M) 

IF ( LTOPT . EQ . 1 ) TOA ( M ) =CT+ULT 
IF(LTOPT.EQ.2)THEN 
CALL ARRIVE 
TOA(M)=CT+ARV 
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END IF 

IF(IP.EQ.1)THEN 

WRITE(7,503)CT,IET,M,IQ(M) ,TOA(M) , (lOH(I) ,1=1 ,4) 

END IF 

IFaP.EQ.2)THEN 

WRITE(7,603)CT,IET,M,IQ(M) ,TOA(M) , lOH(M) 

END IF ■ 

IC = IC + 1 

IF (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2) THEN 
IC = 0 

WRITE (7, 200) 

END IF 
END IF 

NOD(M,ITER) = NOD (M, ITER )+l 
OC(M,ITER) = CC(M,ITER)+OCG+OCI(M) 

IF(KEY.EQ.l)GO TO 20 
IF(KEY.EQ.2)GO TO 92 

END IF 

DON'T ORDER UNTIL PREVIOUS ORDER ARRIVES - — 

IF(M.EQ.IT(M,K) )THEN 
IET=0 

IF(IOH(M) .LT.O) THEN. 

NBO(M,ITER)=-IOH(M) 

BOT(M , ITER)=TOA(M)-CT 

IF(NBO(M,ITER) . EQ • 1 )TBOT(M, ITER)=TBOT (M , ITER)+BOT (M , ITER) 

UBO ( M 7ITER ) =UBO (M , ITER) +BOT (H , ITER) 

IF(IP.E0.1)THEN 

WRITE (7 7504 )CT,IET,N,NBQ(M, ITER) ,BOT(M, ITER) ,UBO(M, ITER) , 

& (lOH(l) ,1=1,4) 

END IF 

IF(IP.EQ.2)THEN 

WRITE (7, 604 )CT,IET,M,NBQ(M, ITER) ,BOT(M, ITER) , UBO (M, ITER) , lOH(M) 
END IF 

IC = IC + 1 

IF (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2) THEN 
IC = 0 

WRITE (7, 200) 

END IF 

ELSE 

UOH(M,ITER)=UOH(M,ITER)+IOH(M)’^DA(M) 

IF(IP.EQ.lkHEN 

WRITE (7, 501 )CT,IET,M,UOH(M, ITER) , (IOH{ I ) , 1=1 ,4) 

END IF 

IF(IP.EQ.2)THEN 

WRITE (7, 601 )CT,IET,M,UOH(M, ITER), lOH(M) 

END IF 

IC = IC + 1 

IF (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2) THEN 
IC = 0 

WRITE(7,200) 

END IF 
END IF 
GO TO 20 
ELSE 

DO 10 J=1,K 

10 IF(IQ(IT{M,J)) .GT.O.OR.IQ(JT(M,J)) .GT.O)GO TO 92 
IET=1 
K=K+1 
IT(M,K)=M 

IF(IOH(M) .LT.0)THEN 

MBQ(M.ITER) = -IOH(M) 

IF(KEY.E0.1)I0(M) = IEO(M)+NBQ(M, ITER) 

IF(KEY.EQ.2)IQ(M) = ME(M)+NBQ(M, ITER) 

IF(LTOPT.EQ.ljTOA(M) = CT + ULT 
IF(LTOPT.EQ.2) THEN 



56 



non 



CALL ARRIVE 
TOA(M)=CT+ARV 
END IF 

BOT(M,ITER) = TOA(K) - CT 

IF(NBQ(M,ITER) .EO.lkBOT(M,ITER)=THOT(M,ITER)+BOT(M,ITER) 
UBO (M , ITER) =UBO (fl , ITER)+BOT (M , ITER) 

IF(IP.EQ.1)THEN- 

WRITE(7,502)CT,IET,M,NBQ(M, ITER), BOT(M, ITER), UBO(M, ITER) , 
& IQ(M),T0A(M),{I0H(I),I=1,4) 

END IF 

IF(IP.EQ.2)THEN 

WRITE(7,602)CT,IET,M,NBQ(M,ITER),BOT(M,ITER) , UBO (M, ITER) , 

& IQ(M) ,TOA{M),IOH{M) 

END IF 

IC = IC + 1 

IF (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2) THEN 
IC = 0 

WRITE (7, 200) 

END IF 
ELSE 

IF(KEY.EQ.1)IQ(M)=IEQ(M) 

IF(KEY.EQ.2)IQ(M)=ME (M)-IOH(M) 

I F ( LTOPT . EQ . 1 J TOA ( M ) =CT+ULT 
IF(LTOPT.EQ.2)THEN 
CALL ARRIVE 
TOA{M)=CT+ARV 

END IF - ' • 

IF(IP.EQ.1)THEN 

WRITE(7,503)CT,IET,M,IQ(M) ,TOA(M) , (IOH( I ) , 1=1 , 4) 

END IF 

IF(IP.EQ.2)THEN 

WRITE(7,603)CT,IET,M,IQ(M) ,TOA(M) ,IOH(M) 

END IF 

IC = IC + 1 

IF (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2) THEN 
IC = 0 

WRITE (7, 200) 

END IF 
END IF 

NOD(M,ITER) = NOD(M,ITER)+l 
OC(M,ITER) = OC(M,ITER)+OCG+OCI(M) 

END IF 

IF(KEY.EQ.l)GO TO 20 

DETERMINE THE JOINTLY REPLENISHABLE ITEM 

91 DO 11 1=1, NIT 

IF (lOH(l) .GT.MR(I).AND.IOH(I) .LE.MC(I) .AND.IQ(I) .EQ.0)THEN 
JT(I,K) = I 
IQ(I) = MEd)-IOH(I) 

I F ( LTOPT . EQ . 1 ) TOA ( I ) =CT+ULT 
IF ( LTOPT. EQ. 2) THEN 
CALL ARRIVE 
TOA(I)=CT+ARV 
END IF 

NODd,ITER)= NOD{I,ITER)+l 
OC(I,ITER) = OC(I ,ITER)+OCI(I) 

JOINT = 1 

IF(IP.EQ.1.0R.IP.EQ.2)THEN 

WRITE(7,506)I,IQ{I),TOA(I) 

END IF 
IC = IC + 1 

IF (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2) THEN 
IC = 0 

WRITE(7,200) 

END IF 
END IF 

11 CONTINUE 

NJTO(ITER) = NJTO(ITER) +JOINT 
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GO TO 20 

92 IF(I0H(M) .LT.O) THEN 

NBQ(M,ITER)=-IOH(M) 

UBO ( M , ITER ) =UBO ( M , ITER ) +B0T ( M , ITER ) 

BOT ( M , ITER ) =TOA (M ) - CT 

IF(NBQ(M,ITER) .EQ. l)TBOT(M, ITER) =TBOT(M, ITER) +BOT(M, ITER) 
IF(IP'.EQ.1)THEN 

WRITE ( 7 , 504 ) CT , lET , M , NBQ (M , ITER) , BOT (M , ITER) , UBO (M , ITER) , 

& (I0H(lhl=l,4) 

END IF 

IF(IP.EQ.2)THEN 

WRITE ( 7 , 604 ) CT , lET , M , NBQ ( M , ITER ) , BOT (M , ITER) , UBO (M , ITER ) , lOH ( M ) 
END IF 

IC = IC + 1 

IF (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2) THEN 
IC = 0 

WRITE (7, 200) 

END IF 

ELSE 

UOH(M,ITER)=UOH(M,ITER)+IOH(M)*DA(M) 

IF(IP.EQ.1)THEN 

WRITE (7, 501 )CT,IET,M,UOH(M, ITER), (IOH( I) ,1=1,4) 

END IF 

IF(IP.EQ.2)THEN 

WRITE (7, 601 )CT,IET,M,UOH(M, ITER), lOH(M) 

END IF 

IC = IC + 1 

IF (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2) THEN 
IC = 0 

WRITE (7, 200) 

END IF 
END IF 
GO TO 20 

SET NEW ON HAND LEVEL FOR THE ITEM : ' lA ' 

93 CT=TA 
IET=2 

IOH(IA)=IOH(IA)+IQ(lA) 

DO 12 1=1, K 

IF(IA.EQ.IT(IA,I))THEM 

IT(IA,I)=0 

IQ(IA)=0 

NTBQ ( I A , ITER ) =NTBQ ( I A , ITER ) +NBQ ( lA , ITER ) 

B0T(IA,ITER)=0. 

NBQ(IA,ITER)=0 
TOA(IA)=10>*10 
END IF 

IF(IA.EQ.JT(IA,I))THEN 

JT(IA,K)=0 

IQ(IA)=0 

NTBQ ( lA , ITER) =NTBQ ( lA , ITER ) +NBQ ( lA , ITER) 

BOT(IA,ITER)=0. 

NBQ(IA,ITER)=0 
TOA(IA)=10>*10 
END IF 

12 CONTINUE 

IF(IP.EQ.1)THEN 

WRITE (7, 500 )CT,IET,IA,(IOH(I) ,1=1,4) 

END IF 

IF(IP.E0.2)THEN 
WRITE ( 7 , 600 )CT , lET , lA , lOH ( lA ) 

END IF 

IC = IC + 1 

IF (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2) THEN 
IC = 0 

WRITE (7, 200) 

END IF 
GO TO 20 
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19 

777 

C 
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— COMPUTE INTER-SUMMARY (ANNUAL) 

IF(lP.NE.4nHEN 
WRITE (7, 200) ITER 
END IF 

DO 15 1=1, KIT 

NTDM(ITER)=NTDM(ITER)+NDM(I,ITER) 

NT0D(ITER)=NT0D(ITER)+N0D(I,ITER) 

OCT(ITER)=CCT(ITER)+OC(I,ITER) 

IF(MSV.EQ.1)THEN 

IF(KB(I,ITER) .GT.0)THEN 

FRC(I,ITER)=KB(I,ITER)*1.0/NOD(I,ITER) 

ELSE 

FRC(I,ITER)=0, 

END IF 

SVL(I,ITER)=(1,-FRC(I,ITER))*100. 

ELSE 

IF(NTBQ(I,ITER) .GT.0)THEN 

FRQ(I , ITER)=NTBQ(I , ITER)*1 . 0/NDM(I , ITER) 

ELSE 

FRQ(I,ITER)=0. 

END IF 

SVL(I,ITER)=(1.-FRQ(I,ITER))*100. 

END IF 
CONTINUE 
IF(IP.NE.4)THEN 
DO 16 1=1, NIT 

WRITE(7,507)I ,NDM(I ,ITER) ,NOD(l ,ITER) ,UOH(I,ITER) ,UBO(I,ITER) , 
&NTBQ(I,ITER) , TBOT( I , ITER) , SVL( I , ITER) , IOH(I ) 

CONTINUE 

WRITE (7, 508 )NTDM( ITER) ,NTOD(ITER) 

WRITE ( 7 , 509 )NJTO ( ITER) 

END IF 

- COMPUTE COSTS 

DO 17 1=1, NIT 

HCOST ( I , ITER ) =UOH ( I , ITER ) *H ( I ) 

SCTT(l,ITER)=UBO(I ,ITER)’'(SCT(I)+H(I)) 

SCTU ( I , ITER ) =NTBQ ( I , ITER ) *SCU ( I ) 

SCOST(I , ITER)=SCTU(I , ITER)+SCTT(I , ITER) 

COST(I , ITER)=HCOST(l , ITER)+SCOST(I , ITER)+OC(I , ITER) 

TCOST( ITER )=TCOST(ITER)+COST(I , ITER) 

IF(IP.NE.4)THEN 

WRITE(7,514) 

DO 18 1=1, NIT 

WRITE(7,515) I, HCOST(I, ITER) ,SCTU(I, ITER) ,SCTT(I, ITER) , 

& OC(I , ITER) ,COST(I , ITER) 

WRITE ( 7 , 5 1 6 ) OCT ( ITER ) , TCOST ( ITER ) 

END IF 

DO 19 1=1, NIT 

KBT(I)=KBT(I)+KB(I,ITER) 

CONTINUE 

-- COMPUTE FINAL-RESULTS (MEAN AND STD. DEVIATIONS) 

WRITE(7,510)NSIM 

DO 21 1=1, NIT 

DO 21 J=1,ITER-1 

NTDMS ( I ) =NTDMS ( I ) +NDM ( I , J ) 

NTODS ( I ) =MTODS ( I ) +NOD ( I , J ) 

TUOH ( I ) =TUOH ( I ) +UOH ( I , J ) 

TUBO ( I ) =TUBO ( I ) +UBO ( I , J ) 

TFRT ( I ) =TFRT ( I ) +TBOT ( I , J ) 

TFRO(l)=TFRO(I)+NTBO(I, J)*l .0 
OCTf (I)=OCTT(I)+OC(I, J) 

THCOST ( I ) =THCOST ( I ) +HCOST ( I , J ) 

TSCTT ( I ) =TSCTT ( I ) +SCTT ( I , J ) 
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21 TSCTU(I)=TSCTU(I)-!-SCTU(I,J) 

DC 22 

AMDK ( I ) =NTDKS ( I ) /THAX 
ANOD ( I >=NTODS ( I ) /TAkK 
AUOH ( I ) =TUOH ( I ) /TMAX 
AUBO ( I ) =TU50 ( I ) /TMAX 
ABOT ( I ) =TFRT ( I ) / TMAX 
ABOQ ( I ) =T FRO ( I ) / TMAX 
ANB(I) =KBT(I)^1 .0/NTODS(I) 

AFRQ ( I )=TFRQ ( I ) /NTDMS ( I ) 

IF(MSV.EQ.1T ASVL(I)=(1.-ANB(I))*100. 

IF(MSV.E0.2) ASVL(I)=(1.-AFRQ(I))*100. 

AOC(I)=OCTT(I)/TMAX 
AHCOST( I )-THCOST(I ) /TMAX 
AS CTU ( I ) =TS CTU ( I ) / TMAX 
AS CTT ( I ) =TS CTT ( I ) / TMAX 

ACOST ( I ) =AHCOST ( I ) +ASCTU ( I ) +ASCTT ( I ) +AOC ( I ) 

22 TVC=TVC+ACOST(I) 

DO 23 1=1, ITER-1 

23 NTJTO=NTJTO+NJTO(I) 

IF(NTJTO.GT.O)ANJ = NTJTO/TMAX 
DO 24 1=1, ITER-1 

24 STV=STV+(TC0ST(I)-TVC)**2 
STDVC=SORT(STV/ (ITER-2) ) 

DO 25 1=1, NIT 

25 WRITE(7,511)I,ANDM(I) ,ANOD(I) ,AUOH(I) ,AUB0(I) ,AB0Q(I),AB0T(I) , 

& ASVL(I) 

DO 26 1=1, NIT 
TANDM=TANDM+ANDM( I ) 

TANOD=TANOD+ANOD ( I ) 

TAOC=TAOC+AOC(I) 

26 CONTINUE 

WRITE (7 , 512)TANDM,TANOD 
WRITE(7,513 ANJ 
WRITE(7,517) 

DO 27 1=1, NIT 

27 WRITE(7,515) I , AHCOST(I ) , ASCTU( I ) , ASCTT (I ) , AOC (I ) , ACOST ( I ) 

WRITE(7 ,516)TAOC,TVC 

WRITE(7,518)STDVC 

STOP 

97 FORMAT(2X,' ///,' * INPUT DATA FOR COSTS :',//, ' ITEM 2X UNIT 

& PRICE' ,2X, 'H. RATE' ,2X, 'H. COST' ,3X, ' ORDER COST ' , 9X , ' B/0 COST', 
&/,36X, 'GROUP' ,6X, ' IND . ' ,4X, 'QTY. ' ,5X, 'TIME' ) 

98 F0RMAT(2X,I2,4X,7(F6.2,3X)) 

100 F0RMAT(F7.3,2X,5(I3,2X) ,3(F6.2,2X)) 

101 F0RMAT(I3,3X,F6.2,3X,I2,3X,2(I1,3X)) 

198 F0RMAT(//,5X, '** INPUT DATA : ' , / , 3X, / , 12X , ' * USED DATA FILE NO. =' 
&,I3,/,12X, '* ORDER POLICY OPTION = ' , 13 , / , 12X , ' * SIMIMULATION LIMIT 
&(YRS) =',I3 , /,12X,'^ START LEVEL OPTION = ' , 13 , / , 12X , ' * PRINTOUT 
&OPTION =' ,I3,/,12X, LEAD TIME TYPE = ' , 13 , / , 12X , /// , 7X , ' ORDER 
&POLICY OPTIONS' ,//,10X, '1 : INDIVIDUAL ORDER POLICY ',/, lOX ,' 2 : JO 
&INT ORDER POLICY' ,//,7X, '* START LEVEL OPTIONS ',//, lOX, ' 1 : REGULA 
&R (IND. POLICY W/ E.O.Q ; JOINT POLICY W/ ORDER-UP QTY. ) ' ,/ , lOX, ' 
&2 : RANDOM (REORDER POINT + 1; ORDER-UP POINT) ',//, 7X, ' * PRINTOUT 
&OPTIONS' ,//,10X, '1 : PRINT ALL ITEM ' , / , lOX, ' 2 : PRINT ONLY EVENT 
&ITEM' ,/,10X, '3 : PRINT ONLY SUMM.ARY' , / , lOX , ' 4 : PRINT ONLY FINAL R 
&ESULT' ,//,7X, '* LEAD TIME TYPES ',//,10X,'l : COSTANT LEADTIME',/, 
&10X,'2 : RANDOM LE.ADTIME (UNIFORM W/ MEAN) ' , // , 6X , ' * SERVICE LEVEL 
& : P',I1, ' = ' ,I3,1X, '%' ,5X, .AVG. LEAD TIME : ',F6.3,' YRS .','(' , 
&I2, ' MON.) ' ,///,3X, 'ITEM' ,3X, 'YR DEMAND' ,4X, ' EOQ ' ,6X, 'MUST BUY' , 
&3X,'CAN BUY' ,3X, 'ORD. UP QTY' ,3X, ' START LEVEL ', 3X ,' EXP . LT.') 

199 F0RM.AT(3X,I2,6X,F7.3,5X,4(I3,3X) ,3X,I3,7X,F7.3) 

200 F0RMAT('l' ,//,35X, SIMULATION RESULTS **',//, 20X, ' (EVENT TYPE 

& 0 : DEMAND OCCURED 1 : ORDER PLACED 2 : ORDER OTY ARRIVED)' 

&,//,5X, 'EVENT' ,4X, 'ITEM S' ,2X, 'BACK ORDER' ,3X, 'UNIT YRS(QTY) ' ,4X, ' 
&0N ORDER' ,10X, 'ON HAND' ,9X, 'JOINT REPLENISHMENT' ,/,lX, 'TIME' ,2X, 
Se'^TYPE' ,10X, 'QTY. ' ,2X, 'TIME' ,3X, 'ON HAND' ,2X, 'B/0. S' ,3X, 'QTY. ' ,2X, 
&'AR/T. ' ,4X, '1' ,4X, '2' ,4X, '3' ,4X, '4' ,2X, 'ITEM#' ,2X, 'OD. QTY. ' ,2X, '0 
&D.AR/T. ' ) 
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300 


FORMAT! ' 


1 ' , // 




&'AT THE 


YEAR 




&,10X, ' UNIT y: 




&'ON HAND 


' ,3X, 


500 


FORMAT! ' 


' ,F6 


500 


?CR1LAT( ' 


' ,F6 


501 


FORMAT! ' 


' ,F6 


601 


FORMAT! ' 


' ,F6 


502 


format! ' 


' ,F6 




& F6.3,1X 


,4(14 


602 


FORMAT! ' 


' ,F6 



:0X,'*** INTER - SUMMARY (ANNUAL) ***',//, 20X, 
,2X,I3,//,3X, 'ITEM #',3X,'# DEMANDS 3X, ' # ORDERS' 
: ' ,7X, 'BACKORDER' ,8X, 'SERVICE' ,6X, 'F/ OH',/,40X, 
I/O ' ,2X, 'QTY. ' ,3X, 'TIME' ,5X, 'LEVEL(%) ' ) 



503 

603 

504 

604 

506 

507 

508 

509 

510 



511 

512 

513 

514 



515 

516 

517 



518 



& F6.3,7X,I4) 



FORMAT! 

FORMAT! 

FORMAT! 

FORMAT! 



FORMAT 
&' YEARS 



,F6.4,2X,I2,3X,I2,37X,I3,1X,F6.3,1X,4(I4,1X)) 

' ' ,F6.4,2K,I2,3X,I2,37X,I3,1X,F6.3,7X,I4) 

' ' ,F6.4,2X,I2,3X,I2,5X,I3,3X,F6.3,11X,F6.2,14X,4(I4,1X)) 

' ' ,F6.3,2X,I2,3X,I2,5X,I3,3X,F6.3,11K,F6.2,20X,I4) 
FORMAT(86X,I2,6X,I3,5X,F6.3) 
FORMAT(5X,I2,9X,I4,7X,I3,9X,F9.4,2X,F6.3,2X,I3,3X,F6.3,5X,F6.2, 
9X,I4) 

FORMAT (IX, 'TOTAL = ' , 7X , 15 , 5X , 15 ) 

FORMATS IX, 'TOTAL NUMBER OF JOINT REPLENISHMENT =',I5) 

'1 ' ,///,20X, FINAL RESULTS (AVG. ) lOX, ' DURING' , 13 , 2X, 
*** ' ,//,3X, 'ITEM #',3X,'# DEMANDS' ,3X, ORDERS ', 9X, ' AVG 
Sc. UNIT YEAR' ,4X, 'AVG. B/0 ', 5X, ' SERVICE ' ,/,40X, 

&'ON HAND' ,3X, 'B/0. S' ,4X, 'QTY. ' ,2X, 'TIME' ,3X, 'LEVEL' 

FORMAT ( 5X , 12 , 6X, F6 . 2 , 5X, F6 . 2 , lOX, F7 . 3 , IX , F7 . 3 , IX , F6 . 2 , 2X , F7 . 3 , 
&2X,F6.'2) 

F0RMAT(1X, 'TOTAL AVG.=' , F7 .2 , 5X, F6 .2) 

FORMAT! IX, 'AVG. NUMBER OF JOINT REPLENISHMENT =',F6.3) 
FORMAT(1X,///,10X, ' « COMPUTED ANNUAL COSTS ($)',//, 2X, 

&'ITEM' ,3X, 'HOLDING' ,9X, 'BACKORDER' ,9X, 'ORDER' ,6X, 'TOTAL' ,/,22X, 'QT 
&Y. ' ,7X, 'TIME' ) 

FORMAT(2X,I2,5X,5(F7.3,4X)) 

FORMAT ( IX, /,33X, 'TOTAL = ' ,2(F8.3 ,3X) ) 

FORMAT! IX, /// ,10X, ' COMPUTED AVG. ANNUAL COSTS ($)',//,2X, 
S'lTEM' ,3X, 'HOLDING' ,9X, 'BACKORDER' ,9X, 'ORDER' ,6X, 'TOTAL' ,/,22X, 'QT 
&Y. ' ,7X, 'TIME' ) 

FORMAT ( IX, /,33X, 'STAND. DEVIATION = ',F8.3) 

END 



C 



C 

c 

c 



c 

c 



SUBROUTINE DEMAND 

DIMENSION ISEED(IO) ,DM(10) ,DA(100) ,MR(10) ,ME( 10 ) , IEQ( 10 ) 
COMMON ICASE , ISEED , DM , DA , ARV , ULT , I , ISTP , MR , ME , KEY , lEQ 
CALL LRND ( ISEED ( ICASE ), U, 1 , 1 , 0 ) 

DA(ICASE) = -(l./DM(ICASE))’''ALOG(U) 

RETURN 

END 



SUBROUTINE ARRIVE 

DIMENSION ISEED(IO) ,DM(10),DA(100),MR(10) ,ME(10) ,IEQ(10) 
COMMON ICASE , ISEED , DM , DA , ARV , ULT , I , ISTP , MR , ME , KEY , lEQ 
CALL LRND(ISEED(I) ,U, 1,1,0) 

ARV= 2.*ULT*U 

RETURN 

END 



C 



SUBROUTINE START 



DIMENSION ISEED(10),DM(10) ,DA(100),MR(10) ,ME(10) ,IE( 
COMMON ICASE , ISEED , DM , DA , ARV , ULT , I , ISTP , MR , ME , KEY , II 
IF(KEY.EQ.l) THEN 

CALL LRND(ISEED(I),U, 1,1,0) 

IS= (IEQ(i)-(MR(l)+l))*U 
ISTP= IS+MR(I)+1 
RETURN 



'( 10 ) 

Q 
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c 

c 



iO 

31 

32 



33 



34 

35 

36 



37 

39 



40 



END IF 

IF (KEY. EQ. 2^ THEN 

CALL LRND(I5EED(I),U, 1,1,0) 
IS= (HE(I)-(MR(I)+1))*U 
:STP= IS+I1R(I)+1 
RETURN 
END IF 
END 



SUBROUT INE OPT ION ( ID , RATE 1 , RATE 2 , KE Y , LTOPT , NS IN , I STA , IP ) 

WRITE(6,35) 

READ(5,34) ID 
WRITE (6, 39) 

READ ( ^ , 40 ) RATE 1 , RATE2 
WRITE(6,31) 

READ (5, 34) KEY 
WRITE (6, 37) 

READ (5, 34) LTOPT 
WRITE (6, 36) 

READ (5, 30) NSIM 
WRITE?6,32) 

READ (5, 34) ISTA 
WRITE (6, 33) 

READ(5,34) IP 
RETURN 
FORITAT(I2) 

FORMAT('l',' * ENTER OPTION FOR ORDER POLICY (1 - 2) ! ',/,5X, ' 
&1 = INDIVIDUAL ORDER POLICY',/, 5X,'2 = JOINT ORDER POLICY ') 
rORMATCl',' * ENTER OPTION FOR START LEVEL (1-2) ! ',/,5X,'l 
&= REGULAR (IND. ORDER POLICY : E.O.Q / JOINT ORDER POLICY : ORDER U 
&P QTY.)’ ,/,5X, '2 = RANDOM BETWEEN (MUST ORDER LEVEL+1) AND (ORDER 
&UP LEVEL) ' ) 

FORMATCl',' * ENTER OPTION FOR PRINTOUT (1 - 4) ! ' , / , 5X , ' 1 = PRI 
&NT ALL ITEM' ,/,5X, ' 2 = PRINT ONLY EVENT ITEM' ,/ , 5X, ' 3 = PRINT ONL 
&Y SUMMARY ',/,5X,'4 = PRINT ONLY FINAL RESULT') 

FORMAT (II) 

FORMAT('l',' * ENTER OPTION FOR DATA FILE (1-4) ! ',/, 

&5X, '1 = DATAl' ,3X, '2 = DATA2' ,3X, '3 = DATA3 ' ,3X, '4 = DATA4' ) 
FORMATCl',' * ENTER SIMULATION LIMIT YEARS !',/,5X,'01 = SIMU 
&LATION LIMIT IS 1 YEAR ' , / , 5X, ' 02 = SIMULATION LIMIT IS 2 YEARS' 

& ,/,7X,' ',/,5X,'99 = SIMULATION LIMIT IS 99 YEARS ') 

FORMATCl',' * ENTER THE TYPE OF LEADTIME(1 - 2) !',/,5X,'l = CON 
&STANT LEADTIME' ,/,5X, '2 = RANDOM (UNIFORM WITH MEAN)') 
FORMATCl',' * ENTER RATE OF STOCKOUT COST TO UNIT PRICE !',/,5X 
&,'1.25 0.91 = STOCKOUT COST PER UNIT IS 125% OF UNIT PRICE ',/, 17X , 
&' STOCKOUT COST PER. TIME IS 91% OF UNIT PRICE') 

F0RMAT(2(F4.2,1X)) 

END 



1. DATA FORMAT ( P2 = 0.99, for items in group 1 ) 



004 050.00 

290.000 025 

041.000 003 

077.000 007 

122.000 010 



99 


2 


1 


025 


112 


177 


-1 


035 


095 


005 


048 


093 


007 


077 


150 



159 010.00 

143 010.00 

109 010.00 

178 010.00 



006.90 000.20 

001.20 000.20 

003.90 000.20 

002.30 000.20 
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2. EXEC PROGRAM 



&TRACE OFF 
&FN = SIM 

&FNO = &CONCAT OF ocFN OUTPUT 

&TYPE Do you need to compile your program ? (Y) 

&READ VAR &R_COMPILE 

&IF &R_COMPILE NE Y &GOTO -RUN 

-H FORTVS &FN 

&IF &RC EQ 0 &GOTO -RUN 

&TYPE Your program did not compile; check for errors. 

&TYPE Do you wish to view the program LISTING file? (Y) 
&READ VAR &RSP1 

&IF &RSP1 EQ Y BROWSE &FN LISTING A 

&TYPE Do you wish to XEDIT the program file? (Y) 

&READ VAR &RESP1 

&IF &RESP1 NE Y ScEXIT 1 

&COMMAND XEDIT &FN FORTRAN A 

&TYPE Do yoa wish to run the program again? (Y) 

&READ VAR &RESP2 

&IF &RESP2 EQ Y &GOTO -H 

&EXIT 1 

-RUN 

FILEDEF 01 DISK &FN DATAl A1 

FILEDEF 02 DISK &FN DATA2 A1 

FILEDEF 03 DISK &FN DATA3 A1 

FILEDEF 04 DISK &FN DATA4 A1 

FILEDEF 07 DISK &FN OUTPUT A1 (LRECL 133 

LOAD &FN (START 

&IF &RC EQ 0 &SKIP 9 

&TYPE Your program did not run correctly; check for errors. 
&TYPE Do you wish to XEDIT the program rile? (Y) 

&READ VAR &RESP3 

&IF &RESP3 NE Y &EXIT 2 

&COMMAND XEDIT &FN FORTRAN A 

&TYPE Do you wish to run the program again? (Y) 

&READ VAR &RESP4 

&IF &RESP4 EQ Y &G0T0 -H 

&EXIT 2 

&TYPE YOUR OUTPUT IS IN THE FILE &FN OUTPUT A 
&TYPE Do you wish to BROWSE your output? (Y) 

&READ VAR &RESP 

&IF &RESP EQ Y &C0MMAND BROWSE &FN OUTPUT A 
&TYPE Print your output file? (Y) 

&READ VAR &RESP7 

&IF &RESP7 EQ Y &C0MMAND PRINT &FN OUTPUT A 
-REDO 

&TYPE Do vou wish to XEDIT the program file? (Y/N) 

&READ VAR*&RESP5 

&IF &RESP5 EQ Y XEDIT &FN FORTRAN A 

&TYPE Do you wish to run the program again? (Y) 

&READ VAR &RESP6 

&RESP56 = &C0NCAT OF &RESP5 &RESP6 
&IF &RESP56 EQ YY &G0T0 -H 
&IF &RESP6 EQ Y &G0T0 -RUN 
&EXIT 



63 



3. RUIJ EXAMPLE 



sim 

Do you need to compile your program ? (Y) 

?S FORTRAN COMPILER ENTERED. 12:11:11 

**MAIN** END OF COMPILATION 1 ****** 

**DEMAND** END OF COMPILATION 2 ****** 

**ARRIVE** END OF COMPILATION 3 ****** 

**START** END OF COMPILATION 4 ****** 

**OPTION** END OF COMPILATION 5 ****** 

VS FORTRAN COMPILER EXITED. 12:11:23 



EXECUTION BEGINS... 

* ENTER OPTION FOR DATA FILE (1-4) ! 

1 = DATAl 2 = DATA2 3 = DATA3 4 = DATA4 

1 



* Enter rate of stockout cost to unit price i 

1.25 0.91 = STOCKOUT COST PER UNIT IS 125% OF UNIT PRICE 
STOCKOUT COST PER TIME IS 91% OF UNIT PRICE 



0.40 0.01 



* ENTER OPTION FOR ORDER POLICY (1-2) ! 

1 = INDIVIDUAL ORDER POLICY 

2 = JOINT ORDER POLICY 



* ENTER THE TYPE OF LEADTIME (1 - 2) ! 

1 = CONSTANT LEADTIME 

2 = RANDOM (UNIFOPJ-I WITH MEAN) 



* ENTER SIMULATION LIMIT YEARS ! 

01 = SIMULATION LIMIT IS 1 YEAR 

02 = SIMULATION LIMIT IS 2 YEARS 

99 ^"simulation LIMIT IS 99 YEARS 
50 



* ENTER OPTION FOR START LEVEL (1-2) ! 

1 = REGULAR (IND. ORDER POLICY: E.0.0 / JOINT ORDER POLICY: OD.UP QTY.) 

2 = RANDOM BETWEEN (MUST ORDER LEVEL+1) AND (ORDER UP LEVEL) 



* ENTER OPTION FOR PRINTOUT (1 - 4) ! 

1 = PRINT ALL ITEM 

2 = PRINT ONLY EVENT ITEM 

3 = PRINT ONLY SUMMARY 

4 = PRINT ONLY FINAL RESULT 

4 

YOUR OUTPUT IS IN THE FILE SIM OUTPUT A 
Do you wish to BROWSE your output? (Y) 

^rint your output file? (Y) 
n 

Do you wish to XEDIT the program file? (Y/N) 
n 

Do you wish to run the program again? (Y) 
n 

R; T=4.63/5.82 12:11:54 
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4. OUTPUT 1 (individual order policy) 



** INPUT DATA : 

* USED DATA FILE NO. = 1 

* ORDER POLICY OPTION = 1 

* SIMIMULATION LIMIT(YRS) = 50 

* START LEVEL OPTION = 2 

* PRINTOUT OPTION = 4 

* LEAD TIME TYPE = 1 



* ORDER POLICY OPTIONS 

1 : INDIVIDUAL ORDER POLICY 

2 ; JOINT ORDER POLICY 

* START LEVEL OPTIONS 

1 ; REGULAR (IND. POLICY - E.O.Q ; JOINT POLICY - OD.UP QTY.) 

2 : RANDOM (REORDER POINT + 1; ORDER-UP POINT) 

* PRINTOUT OPTIONS 

1 : PRINT ALL ITEM 

2 : PRINT ONLY EVENT ITEM 

3 : PRINT ONLY SUMMARY 

4 : PRINT ONLY FINAL RESULT 

* LEAD TIME TYPES 



1 : COSTANT LEADTIME 

2 : RANDOM LEADTIME (UNIFORM W/ MEAN) 



* SERVICE LEVEL: 


P2= 99 


% * 


AVG. LEAD 


TIME: 0.083 YRS . ( 


1 MON.) 


ITEM 


YR DEMAND 


EOQ 


MUST BUY 


CAN BUY 


ORD . UP START 


EXP. LT. 


1 


290.000 


159 


25 


112 


177 65 


0.083 


2 


41.000 


143 


3 


35 


95 45 


0.083 


3 


77.000 


109 


7 


48 


93 38 


0.083 


4 


122.000 


178 


10 


77 


150 60 


0.083 



* INPUT 


DATA FOR 1 


COSTS: 








ITEM 


UNIT PRICE 


H.RATE H.COSl 


ORDER COST B/0 COST 










GROUP 


IND. QTY. 


TIME 


1 




6.90 


0.20 1.38 


50.00 


10.00 2.76 


0.07 


2 




1.20 


0.20 0.24 


50.00 


10.00 0.48 


0.01 


3 




3.90 


0.20 0.78 


50.00 


10.00 1.56 


0.04 


4 




2.30 


0.20 0.46 


50.00 


10.00 0.92 


0.02 






** FINAL RESULTS (AVG. 


) : 


DURING 50 YEARS 


*** 


ITEM 


It 


# DEMANDS 


# ORDERS 


AVG. UNIT YEAR AVG. B/0 


SERVICE 








ON HAND 


B/O.S QTY. TIME 


LEVEL (%) 


1 




292.24 


1.84 


79.151 


0.025 2.92 0.009 


99.00 


2 




41.58 


0.30 


69.217 


0.008 0.34 0.006 


99.18 


3 




77.00 


0.72 


54.064 


0.004 0.46 0.003 


99.40 


4 




120.50 


0.68 


88.760 


0.014 0.84 0.006 


99.30 


TOT AVG. 


; 531.32 


3.54 








AVG. 


NO. 


OF JOINT 


REP; 0.000 
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ITEM 



* COMPUTED AVG. AMNUAL COSTS ($) 



HOLDING 


BACKORDER 


ORDER 


TOTAL 




QTY. 


TIME 






109.228 


8.059 


0.036 


110.400 


227.723 


16.612 


0.163 


0.002 


18.000 


34.777 


42.170 


0.718 


0.003 


43.200 


86.090 


40.830 


0.773 


0.007 


40.800 


82.409 






TOTAL 


= 212.400 


430.999 






STAND. 


DEVIATION = 


51.942 
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5. OUTPUT 2 (joint order policy) 



** INPUT DATA : 

* USED DATA FILE NO. = 1 

* ORDER POLICY OPTION = 2 

* SIMIMULATION LIMIT(YRS) = 50 

* START LEVEL OPTION = 2 

* PRINTOUT OPTION = 4 

* LEAD TIME TYPE = 1 



* ORDER POLICY OPTIONS 

1 : INDIVIDUAL ORDER POLICY 

2 : JOINT ORDER POLICY 

* START LEVEL OPTIONS 

1 : REGULAR (IMD. POLICY - E.O.Q ; JOINT POLICY - OD. UP QTY.) 

2 ; RANDOM (REORDER POINT + 1; ORDER-UP POINT) 

* PRINTOUT OPTIONS 

1 : PRINT ALL ITEM 

2 : PRINT ONLY EVENT ITEM 

3 : PRINT ONLY SUMMARY 

4 : PRINT ONLY FINAL RESULT 

* LEAD TIME TYPES 



1 : COSTANT LEADTIME 

2 : RANDOM LEADTIME (UNIFORM W/ MEAN) 



SERVICE LEVEL: 


P2= 99 


% * AVG. 


LEAD TIME 


: 0.083 


YRS.( 1 


MON. ) 


ITEM 


YR DEMAND 


EOQ 


MUST BUY 


CAN BUY 


ORD. UP 


START 


EXP. LT 


1 


290.000 


159 


25 


112 


177 


71 


0.083 


2 


41.000 


143 


-1 


35 


95 


23 


0.083 


3 


77.000 


109 


5 


48 


93 


32 


0.083 


4 


122.000 


178 


7 


77 


150 


50 


0.083 



* INPUT DATA FOR COSTS: 



ITEM 


UNIT PRICE 


H.RATE H.COSI 


■ ORDER COST 


B/0 COST 








GROUP 


IND. 


QTY. 


TIME 


1 


6.90 


0.20 1.38 


50.00 


10.00 


2.76 


0.07 


2 


1.20 


0.20 0.24 


50.00 


10.00 


0.48 


0.01 


3 


3.90 


0.20 0.78 


50.00 


10.00 


1.56 


0.04 


4 


2.30 


0.20 0.46 


50.00 


10.00 


0.92 


0.02 




** FINAL RESULTS (AVG. 


) : 


DURING 50 YEARS 




ITEM 


# # DEMANDS 


# ORDERS 


AVG. UNIT YEAR 


AVG. B/0 


SERVICE 






ON HAND 


B/O.S 


QTY. TIME 


LEVEL (% 


1 


292.24 


1.94 


75.843 


0.021 


2.56 0.007 


99.12 


2 


41.58 


0.60 


56.035 


0.000 


0.00 0.000 


100.00 


3 


77.00 


1.13 


51.541 


0.001 


0.04 0.001 


99.95 


4 


120.50 


1.06 


30.055 


0.003 


0.30 0.001 


99.75 



TOT AVG.; 531.32 4.78 

AVG. NO. OF JOINT REP; 1.620 



67 



* COMPUTED AVG. ANNUAL COSTS ($) 



ITEM 


HOLDING 


BACKORDER 


ORDER 


TOTAL 


1 


104.663 


QTY. 

7.066 


TIME 

0.030 


113.400 


225.159 


2 


13.448 


0.000 


0.000 


6.000 


19.448 


3 


40.202 


0.062 


0.001 


12.800 


53.065 


4 


36.825 


0.276 


0.001 


13.600 


50.703 








TOTAL 
STAND . 


= 145.800 

DEVIATION = 


348.375 

21.055 
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